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I,    i    =   ideal  value 

Yd,   ye  -  system  error 

Q,   q    ss   Indirectly  controlled  variable 
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*  Above  are  standard  symbols  taken  from  reference  1. 
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Chapter  1 
THE  PROBLEM 
1.1   General* 

Feedback  Control  System  synthesis  means  the 
determination  of  system  and  component  specifications 
to  meet  the  requirements  of  a  specified  job.   The  first 
phase  of  such  a  procedure  involves  the  selection  or 
design  of  a  power  element  adequate  to  drive  the  load.   This 
paper,  however,  is  not  concerned  with  the  selection  of 
components,  and  the  power  element  is  stipulated. 

The  type  of  Feedback  Control  System  considered  is  a 
servomechanism.   By  definition  here,  a  servomechanism  is  a 
particular  type  of  feedback  control  system  in  which  the 
controlled  variable  is  a  mechanical  position.   The  out- 
put is  the  mechanical  position  of  one  object  relative  to 
another. 

The  words  "synthesis"  and  "analysis"  are  frequently 
used  interchangeably  and  loosely.   Here,  we  will  use 
J.  R.  Burnett's  "  definitions 

(1)  The  synthesis  problem.   Given  the  input  to  a 
system  and  the  required  output,  determine  the  transfer 
function  of  the  system. 

(2)  The  analysis  problem.   Given  the  input  to  a 
system  and  the  transfer  function  of  the  system,  find  the 

*  numerical  superscripts  refer  to  Bibliography 


output  of  the  system. 

The  type  of  system  under  consideration  is  lumped 
parameter,  finite,  time- invariant,  and  linear* 

Stability  is  a  problem  to  be  reckoned  with*  We  can 
say  that  stability  is  the  primary  consideration  in  all 
control  systems.   In  that  sense,  it  can  be  considered  as  a 
basic  specification,  always  implied.   However,  in  another 
sense,  it  is  of  only  secondary  consideration,  because  it  is 
virtually  always  possible  to  render  a  system  stable  by  some 
form  of  compensation.   For  linear  systems,  stability  is  a 
function  of  the  system  alone  and  is  not  dependent  upon  the 
input  to  the  system.   A  system  whose  response  will 
eventually  become  arbitrarily  small,  once  the  input  is 
removed  is  "stable".   Stated  another  way,  a  system  is 
stable,  if  the  impulse  response  of  the  system  approaches  a 
constant  value  and  remains  constant  for  large  values  of 
time  after  the  impulse  has  occurred.   If  the  closed-loop 
system  function  is  stated  as 

R(*0     s^bm.,S% D(SJ 

then  D(s)  is  a  Hurwitz  polynomial,  all  of  whose  roots  must 
have  a  negative,  non-zero,  real  part,  for  stability. 

In  the  preparation  of  this  paper  we  have  researched 
a  great  many  writings.  We  find  all  of  them  in  agreement 
on  one  thing:   In  the  design  of  Feedback  Control  Systems, 
arriving  at  an  accurate  and  complete  statement  of  the 


problem  is  the  most  difficult  phase.   The  difficulty  of 
this  phase  is  aggravated  by  the  fact  that  Feedback  Control 
Systems  have  permeated  all  walks  of  our  modern  life.   The 
science  of  feedback  control  knows  no  bounds  as  it  cuts 
across  all  historic  lines  of  endeavor.   The  science  has 
become  so  broad  with  so  many  people  affected  by  it,  that 
a  statement  of  a  particular  design  problem  has  little 
universal  application.  Each  problem  statement  seems  to  be 
unique  unto  itself. 

The  vastness  of  this  new  science  is  evidenced  by  the 
fact  that  many  engineers  no  longer  call  themselves 
"Electrical"  or  "Mechanical"  or  "Aeronautical"  Engineers, 
but  rather  "Control"  Engineers.   Many  books  are  being 
written  and  published  as  "Control  Engineering"  books. 
There  is  a  great  deal  of  discussion  encouraging  the  divorce 
of  Feedback  Control  Systems  from  all  other  historical  lines 
of  engineerings  to  let  it  grow  as  an  independent  science. 
1.2  Flow  Diagram  of  Design 

A  typical  design  flow  diagram  might  be  as  follows? 
(See,  for  example,  writings  of  R.  B.  Wilcox  ). 
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Note  that  "feedback"  continually  exists  through  all  the 


design  phases.   This  is  what  makes  control  system  design 
so  fascinating,  but  at  the  same  time  so  maddening. 

The  problem  statement  may  be  given  in  any  one  of 
three  forms; 

(1)  a  literal  description  of  what  the  system  is 
supposed  to  do,  or 

(2)  a  numerical  specification  of  the  system  per- 
formance ,  or 

(3)  a  graphical  representation  of  the  system. 
From  the  designees  viewpoint,  numerical  specification 
is  preferred  to  eliminate  ambiguity,  but  the  price  is 
loss  of  flexibility  in  meeting  changing  design  considera- 
tions.  If  no  graphical  representation  is  provided  the 
designer  is  immediately  required  to  prepare  a  functional 
block  diagram  of  some  kind  so  that  he  may  proceed  with  a 
component  design  analysis.   Completion  of  this  phase  leads 
to  the  determination  of  suitable  transfer  functions  for  the 
components.   Component  transfer  functions  may  be  determined 
by  experimental  or  analytical  means.   There  are  two  basic 
experimental  methods;  determine  the  transfer  function  from: 

(1)  the  component  frequency  response 

(2)  the  component  transient  response. 
Analytical  methods  require  the  use  of  basic  mathematical 
and  physical  principles  (for  example,  Newtonfs  Laws  of 
Motion) • 

In  any  event,  a  system  block  diagram  may  then  be 
constructed,  from  which  the  "system"  (or  what  we  prefer  to 
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here  call  "plant")  transfer  function  is  determined  in 
either  algebraic  or  graphical  form.   A  pitfall  comes 
to  light  here,  because  at  this  juncture,  there  is  a 
tendency  to  forget  the  limitations  imposed  on  the 
assumed  mathematical  model.   These  limitations  might 
result  from  approximation,  initial  conditions,  noise, 
drift  or  some  other  special  input. 

System  analysis  might  then  proceed  using  time  domain 
or  frequency  domain  methods  with  or  without  the  aid  of 
analog  or  digital  computers.   It  is  this  analysis  which 
then  forms  a  basis  for  any  synthesis  or  modification 
technique. 
1.5  Problem  Statement 

The  type  or  manner  of  statement  of  the  problem 
requirements  often  determines  the  feasibility  of  any 
particular  approach.   Although  ideally,  the  requirements 
are  documented  in  the  form  of  complete  numerical  specifi- 
cations, more  often  than  not,  the  only  information  is  the 
general  requirement  that  some  given  operation  must  be 
automatically  controlled.   The  designer,  therefore, 
must  conduct  a  study  to  formalize  his  own  detailed  specifi- 
cations. 
1.4  Component  Design  and  Analysis 

This  paper  is  not  concerned  with  this  phase  of 
design.   It  will  here  be  assumed  that  the  components  of 
the  plant  are  specified.  We  are  concerned  with  the 


component  design  of  the  compensator  only  in  the  sense 
of  whether  the  component  is  physically  realizable. 

1.5  Transfer  Functions 

Since  we  here  assume  linear  time -invariant  systems, 
transfer  functions  may  be  manipulated  following  algebraic 
rules;  isolation  between  component  transfer  functions 
is  implied.   It  follows  that  the  concept  of  superposition 
is  valid.   This  means  the  output  of  the  system  may  be 
determined  as  the  result  of  all  the  various  inputs  and 
disturbances,  each  applied  separately.   It  is  also 
possible  to  distinguish  between  dynamic  and  steady- state 
performance  in  each  of  these  cases.   The  error  may  also 
be  analyzed  both  dynamically  and  in  steady  state. 
However,  a  transfer  function  represents  the  assumption  of 
some  mathematical  model,  and  as  such,  the  limitations 
imposed  by  the  assumptions  must  be  kept  constantly  in 
mind.   In  the  mechanics  of  attempting  system  optimization, 
for  example,  the  limitation  of  saturation  is  always 
present.   As  the  signal  approaches  saturation,  the  assump- 
tion of  linearity  becomes  less  and  less  valid. 

1.6  System  Analysis  and  Synthesis 

The  time  and  frequency  analysis  of  the  system 
characteristics  is  an  analytical  and/or  graphical  predic- 
tion that  the  system  will  perform  as  the  dynamic  specifica- 
tions require.   If  it  does  not,  then  synthesis  is 
required  to  determine  necessary  modifications  and 
revisions  to  fulfill  the  response  and  accuracy 
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specifications.   Synthesis  and  Analysis  then  go  hand  in 
hands  the  final  analysis  is  the  one  which  defines  all  the 
characteristics  of  the  completed  design.  We  generally 
refer  to  system  modifications  and  revisions  as  system 
"compensation".   Compensation  will  mean  the  inclusion  of 
that  network  which  is  necessary  in  order  tos   (1)  make 
the  system  stable ,  and  (2)  meet  the  specifications. 
A  study  of  the  literature  reveals  this  phase  to  be  the  one 
where  writers  and  designers  reach  a  parting  of  the  ways. 
They  proceed  down  one  of  two  broad  paths.  The  best 
descriptive  titles  for  these  paths  that  we  have  found  are 
those  used  by  Newton^s 

(1)  Trial  and  Error  Design  Methods,  and  (2)  Analytical 
Design  Methods.   Some  of  the  more  lucid  writings  in  the 
current  literature  on  these  two  philosophies  are  those  of 
Thaler  and  Brown^  for  Trial  and  Error  Design  Methods  and 
those  of  Newton^  for  Analytical  Design  Methods. 
1.7  Objective  of  This  Paper 

In  this  paper,  we  attempt  to  bring  under  one  cover 
a  synopsis  of  the  current  methods  suggested  or  in  use  for 
the  design  of  the  servome onanism.   Our  objective  is  to 
provide  insight  into  that  phase  of  design  surrounding  the 
problem  statement  where  the  designer  makes  his  decision  as 
to  what  approach  to  attempt  first.   To  that  end  an  under- 
standing of  the  capabilities  and  limitations  of  the 
various  techniques  is  required.   The  methods  considered, 
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although  not  all  inclusive ,  are  those  which,  in  the 
opinion  of  the  authors ,   do  have  the  greatest  usefulness 
or  offer  the  most  potential • 

In  Chapter  2  we  explore  more  thoroughly  the  statement 
of  specificationso   Chapters  3  and  4  contain  a  summation 
of  the  basic  procedures  and  mechanics  employed  in  the  above 
two  design  methods*  In  Chapters  5  and  69  we  have  listed 
some  of  the  modifications  and  refinements  that  are  found  in 
the  current  literature.   The  compilation  is  not  all 
inclusive  but  the  techniques  listed  attempt  in  some  way  to 
make  the  problem  of  compensation  easier  within  their 
specified  limitations. 
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CHAPTER  2 
SPECIFICATIONS 
2.1  General 

It  may  often  happen  that  the  specifications  given  to 
the  servo  designer  are  either  incomplete,  incompatible .  or 
incomprehensible*   Specifications  fall  into  two  main 
categories?   (1)  The  specification  of  control- system 
dynamics  and  performance ,  and  (2)  the  general  specifications^ 
such  as  power- supply  variations  and  environmental  conditions 5 
which  influence  the  dynamic  characteristics  and  performance* 
We  mentioned  earlier  that  the  designer  usually  must  conduct 
a  study  to  formalize  his  own  detailed  specifications. 
This  study  may  immediately  reveal  contradictory  specifica- 
tions.  It  may  indicate  that  relaxation  of  a  particular 
requirement  will  greatly  simplify  the  control  system.   It 
may  indicate  that  one  specification  is  so  domineering  that 
satisfaction  of  that  one  is  tantamount  to  overall  success. 

A  servomechanism  is  expected  to  perform  any  or  all  of 
the  following  functions s 

(1)  Bring  about  a  change  in  the  actual  value  of  the 
output  so  that  it  conforms  to  a  desired  value  at  all 
times 5 

(2)  Minimize  the  effect  of  varying  component  per- 
formance on  the  outputs 

(3)  Minimize  the  effect  of  disturbances. 

The  performance  specifications  determine  the  degree  of 
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excellence  with  which  the  servo  must  carry  out  the  above 
functions*   This  means  that  the  performance  specifications 
must  be  given  in  terms  of  the  desired  output  for  a  given 
inputo   It  is j,  therefore 9   necessary  to  explore  the  type  of 
inputs  that  might  arise • 
2.2  Inputs 

The  input  signal  may  be  one  of  the  following  types? 

(a)  Aperiodic 5  noise  free 

(b)  Aperiodic p  with  noise 

(c)  Periodic,  noise  free 

(d)  Periodic ,  with  noise 

(e)  Stochastic 9   noise  free 

(f)  Stochastic j  with  noise 

(Noise  is  regarded  as  any  input- signal  variation  that  is 
not  a  measure  of  the  information  carried  by  the  input.) 

Some  typical  specifications  for  the  six  types  of 
input  signals  listed  are  as  follows  (see  chapter  16  of 
reference  4)s 

(a)  Aperiodic s   noise  free 

(1)  The  system  dynamic  error  shall  neither  exceed 
a  specified  maximum  value s   nor  shall  the 
steady  state  error  exceed  a  specified  maximum 
value • 

(2)  The  integral- square  system  error  shall  not 
exceed  a  specified  value. 
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(3)  The  system  «rror?  in  addition,  shall  not  exceed 
a  specified  amount  in  the  presence  of  a 
specified  disturbance  that  occurs  at  some 
specified  point  in  the  systemc 

(4)  See  paragraph  4«2  for  other  possible  error 
specificationso 

(b)  Aperiodic,  with  noise 

(1)  Given  that  the  first  component  of  the  system 
error  is  that  for  zero  noise  step  inputs  the 
second  component  is  the  value  of  the  output 
from  noise  alone j  then  the  square  root  of  the 
sum  of  the  squares  of  the  two  components  shall 
not  exceed  a  specified  value* 
(Aperiodic  signals  commonly  considered  are  steps,  ramps, 
impulses 5  pulses,  or  an  input  expressed  as  a  power  series 
in  time*) 

(c)  Periodic,  noise  free  (only  fundamental  frequency 
present) 

(1)  The  frequency  response  shall  be  characterized 
by  a  specified  peak  magnitude  ratio  (output/ 
input)  occurring  at  a  specified  frequency*, 

(2)  The  magnitude  ratio  (output/input)  shall  be 
within  a  band  of  some  specified  number  of 
decibels  over  a  specified  frequency  range,  and 
phase  shift  (output/input)  shall  not  exceed  a 
specified  amount  over  this  same  range • 
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(d)  Periodic,  with  noise 

(1)  Error  expressed  in  a  similar  manner  to  that  of 
aperiodic,  with  noise,  above* 

(e)  Stochastic,  noise  free  (see  paragraph  2*4  for 
definition.) 

(1)  System  error  shall  not  exceed  a  specified  rms 
value  when  the  input  autocorrelation  function 
has  a  given  value « 

(f)  Stochastic,  with  noise 

(1)  System  error  shall  not  exceed  a  specified  rms 
value  when  the  input  signal  autocorrelation 
function,  the  input-noise  autocorrelation 
function,,  and  the  signal°to-noise  cross- 
correlation  function  are  given* 
2.5  Disturbances 

It  is  also  necessary  to  specify  performance  in 
response  to  a  given  load  and/or  disturbance  occurring  at 
points  different  from  the  input*   The  load  or  disturbance 
can  also  be  classed  as  aperiodic,  periodic,  or  stochastic* 
Typical  specifications  take  on  the  same  form  as  those  above 
for  inputs.   A  load  specification,  however,  usually 
prescribes  the  amount  of  time  allowed  for  the  output  to 
recover  to  within  a  specified  deviation* 
204   Stochastic  Signals 

A  stochastic  process  is  one  in  which  there  is  an 
element  of  chance*   Sometimes  the  input  to  a  system  is  not 
completely  predictable  and  cannot  be  described  by  a 

13 


mathematical  function  of  time.  A  typical  example  of  a 
stochastic  process  is  a  radar  signal  mixed  with  noise. 
Since  the  value  of  a  stochastic  signal  cannot  be  determined 
with  certainty  at  a  given  instant  of  times  probability 
density  functions  and  other  statistical  characterizations 
such  as  the  average  value D    the  rms  value 0    and  the  correla- 
tion function  are  used  to  describe  the  signal  (see  Appendix 
B)  „  However.,  it  is  necessary  to  think  of  a  stochastic 
signal  as  a  member  of  a  family  of  signals.,  each  generated 
by  an  identical  process©  Such  a  family  of  signals  is  called 
an  ensemble  and  the  statistical  characterization  (such  as 
a  correlation  function)  of  the  stochastic  process  is  related 
to  the  ensemble  rather  than  to  a  particular  member  of  the 
ensemble •   It,,  therefore „  follows  that  the  determination  of 
the  response  of  a  system  to  a  stochastic  Input  does  not 
yield  a  function  of  time s   but  rather  a  statistical 
characterication  of  the  output  signal  ensemble  (see  Chapter 
4). 
2c5  Philosophy  on  Choice  of  Test  Input 

It  may  be  necessary  for  the  designer  to  prescribe  his 
own  inputs  when  analyzing  the  effects  of  noise 9    or  load 
disturbances.,  or  the  effects  of  environmental  conditions 
such  as  temperature j,  humidity 9    corrosion,  etc.  Even  though 
the  time  domain  characteristics  are  frequently  specified  in 
terms  of  the  response  ratio  of  the  output  to  a  step- function 
input t   many  others  may  be  specified  or  iinplied.   Further- 
more P  they  may  be  extremely  complicated — requiring 
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graphical  description  either  in  the  time  domain  or  as  a 
power  density  spectrum.  It  might  be  well  here  to  quote 
the  feelings  of  some  writers  in  the  field.  It  is  to  be 
noted  that  they  are  not  all  in  agreement.  For  example 0 
the  following  quotation  (page  308 ,  reference  5)  succinctly 
expresses  one  viewpoint: 

"The  characteristics  of  a  particular  servo  should  be 
determined  by  the  actual  input,  the  actual  uncontrolled 
disturbances  acting  on  the  system  and  the  actual  output 
requirements*   It  is  clearly  not  sufficient  to  assume  the 
input  to  be  a  step  in  displacement  or  velocity,  nor  is  it 
sufficient  to  require  only  that  the  transient  response  be 
well  damped  and  that  the  velocity  lag  be  small  ••<,.  In 
general,  the  actual  input  and  noise  as  well  as  the  output 
requirements  need  a  statistical  description  »,,,   One  can, 
of  course,  conceive  of  specialized  servo  problems,  in  which 
the  input  is  a  displacement  step  and  in  which  the  require- 
ments are  based  on  the  transient  response;  but  they  are  the 
exception  rather  than  the  rule,"   (words  of  R,  S.  Phillips), 

On  the  other  hand,  quoting  from  the  same  reference 
(page  18) s   "The  performance  of  a  servo  can  also  be 
specified  in  terms  of  its  response  to  a  step  function* 
The  procedure  of  experimentally  and  theoretically  studying 
a  servo  through  its  response  to  a  step- function  input  is 
extremely  useful  and  is  widely  used  for  a  number  of 
reasons.   The  experimental  techniques  used  in  such 
testing  are  simple  and  require  a  minimum  of  instrumentation, 
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The    characteristics   of  any  truly  linear   system  are,    of 
course j,    completely  summarized  by  its   response    to   a   step- 
function  input a      That  is,    if  the    step-function  response 
is  known,    the   response    to   any  other  arbitrary  input   signal, 
can  be   determine do      It  would  be   expected,    therefore,    and  it 
is   true,    that  with  proper  interpretation  the    step- function 
response   is   a  powerful  and  useful  criterion  of  overall 
system  quality*,11      (words   of  I.   A*.    Getting). 

But,   progress  and   advances   in  technology  must  be 
recognized©      Less   than  ten  years   later,   Dr»   T,   C,   Fry,   Bell 
Telephone   Lab,    Inc.      saids    "...•    in  any  actual  guidance   or 
fire<°control  problem,   we   are   not  really  concerned  with  the 
response    of  the    system  to   some   particular,    ideally  defined 
tactical   input •     We   are   concerned  with  its   response    to  the 
whole   gamut   of  possible   tactical   situations,    including  all 
the    possible   variability  in  the   enemy  path  and   all   the 
possible   errors  which  may  be   produced  by  random  perturba- 
tions  in  the    input   data  or  in  the   mechanism  itself • 
Obviously,    this  adds  elements   of  information  theory, 
statistical   theory  or  whatnot   and  greatly  increases   the 
level   of   (essentially  mathematical)    insight  required  for 
effective  work," 

The  writers   seem  to  be   hitting  at  the    same   old  saw  of 
the    tug-of-war  between  abstract  mathematics   and 
engineering   approximations*      To  easily  and  practically 
accomplish  a  design  by  hand,    on  paper,    requires   simple 
inputs.      The  use   of  the    step  input,   for  example,   does   take 
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into  consideration  the  the  whole  frequency  range  •   But., 
as  control  systems  become  more  and  more  sophisticated^  it 
becomes  necessary  to  consider  more  sophisticated  inputs. 
This  fortunately  becomes  feasible  in  the  modern  day  with 
the  use  of  computing  machinery* 
206  Static  Characteristics 

The  static  performance  specifications  describe  the 
steady  state  value  of  the  system  output*   Although  the 
statement  of  them  is  a  simple  matters  their  importance 
lies  in  the  fact  that  they  sometimes  immediately  establish 
the  type 3  system  which  must  be  used*  For  example 9    if  a 
system  is  to  have  no  steady  state  error  in  response  to  a 
step  input 5  it  must  be  at  least  Type  1  (meaning  one  pure 
integration  must  exist  in  the  open-loop  transfer  function) £ 
if  no  steady  state  error  is  to  exist  in  response  to  a  ramp 
inputs  the  system  must  be  at  least  Type  2*      (Two  pure 
integrations  must  exist  in  the  open-loop  transfer  functions*) 

Normally  an  acceptable  following  error  must  be  stated 
for  the  response  to  a  ramp  input*   Some  maximum  allowable 
output  is  usually  stated  for  the  response  to  a  disturbing 
input o 
2*7  Dynamic  Characteristics 

The  desired  dynamic  characteristics  may  be  specified 
in  terms  of  transient  or  frequency  response©   (If  defined 
in  both  domains,,  their  compatibility  must  be  determined*) 
The  three  main  themes  of  engineering  design  might  be 
tabulated  as  follows? 
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SPEED  OF  RESPONSE 

RELATIVE  STABILITY 

Transient 
Domain 

Frequency 
Domain 

Transient     Frequency 
Domain       Domain 

a)  time 
constant 

b)  rise 
time 

c)  settling 
time 

a)  bandwidth 

b)  cutoff 
frequency 

c)  frequency  of 
peak 
overshoot 

a)#  of  first 
overshoot 

b)  number  of 
overshoots 

a)  Peak  value 
of  closed 
loop  output- 
input 
amplitude 
ratio 

b)  gain  and 
phase  margin 
of  open- 
loop 

frequency 
characteris- 
tic 

ACCURACY 

Allowable  error  : 

Ln 

terms 

of 

%   or  per  unit 

of 

controlled 

variable9  stated 

as 

9 
o 

a) 

Maximum 

b) 

Ave  rage 

c) 

rms  value 

— —. 
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2<>8      Required  Specifications 

The   three   main  themes  of  design  specifications   are 
essentially  those    of   (1)    speed  of  response,,    (2)    stability,, 
and    (3)    accuracy a      The   Trial   and  Error  Design  Methods 
specification  requirements   ares 

(1)  Input  Signal 

(2)  Desired  Output 

(3)  Disturbances  and  special  inputs 

(4)  Allowable  error 

(5)  Plant  elements  (fixed) 

(6)  Relative  stability 

The  Analytical  Design  Methods  specification  requirements 
ares 

(1)  Input  signal 

(2)  Desired  Output 

(3)  Disturbances  and  special  inputs 

(4)  Performance  index  and  required  value  of  same 

(5)  Degree  of  Freedom  allowed  in  compensation 
Notice  that  the  specifications  for  the  two  methods 

have  two  basic  differences? 

(a)  The  Analytical  Design  Method  calls  for  a 
performance  index  vice  allowable  error  and  relative 
stability  specifications •   A  performance  index  is  simply 
a  single  number  which  is  used  as  an  indirect  measure  of 
system  performance 0   Its  use  is  an  attempt  to  replace  the 
functional  description  of  the  performance  of  a  system 
through  its  response  parameters  (such  as  peak  overshoot,, 
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rise  time,  etc.)  with  a  numerical  description  that  rates 
the  system  performance  with  a  single  number* 

(b)  The  Analytical  Design  Method  calls  for  the  degree 
of  freedom  allowed  in  compensation  vice  relative  stability 
and  plant  elements  specifications*   This  specification  is 
not  strictly  necessary,  but  practically  it  is.  With 
constraints  imposed  by  degrees  of  freedom,  it  is  possible 
to  categorize  large  portions  of  the  computational 
mathematics  for  all  time,  and  make  use  of  only  the  results 
of  this  categorization  in  a  particular  problem*   More  will 
be  said  about  these  differences  in  the  next  two  chapters* 
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CHAPTER  3 
BASIC  TRIAL  AND  ERROR  DESIGN  METHODS 

3.1  General 

Almost  all  Trial  and  Error  Design  Methods  depend 
upon  the  ability  to  express  the  output  of  a  component  or  a 
servome onanism  with  respect  to  an  input  in  terms  of 
differential  equations o   These  equations  are  almost 
universally  transformed  into  algebraic  equations  by  the 
Laplace  transformo   After  this  transformation  is  made, 
components  can  be  collected  into  one  overall  mathematical 
model  by  well  known,  easily  applied  methods,  and  analysis 
and  synthesis  of  the  system  is  made  directly  in  terms  of 
these  equations© 

Most  methods  utilize  and  depend  upon  the  open-loop 
transfer  function  or  equation,  the  closed-loop  transfer 
function,  or  a  combination  of  the  two  functions o   Test 
inputs  with  simple  Laplace  transforms  are  applied  to  these 
equations,  the  output  is  compared  to  the  input,  and  various 
parameters  are  measured  to  determine  the  acceptability  of 
the  mathematical  model  of  the  servomechanismo   Then  the 
necessary  hardware  is  determined  to  match  the  model,  or 
the  model  is  changed  and  another  comparison  or  analysis  is 
made  o 

Specifications  and  design  parameters  which  must  be 
satisfied  by  servomechanism  adjustment  are  expressed  in 
terms  of  either  transient  response  to  a  nonperiodic  test 
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input  or  frequency  response  to  a  periodic  test  input 
(almost  universally  a  sinusoidal  function)* 

Several  design  methods  work  with  the  transfer 
functions  directly  in  the  frequency  domain,,  taking 
advantage  of  the  fact  that  the  analytical  or  graphical 
results  will  be  directly  related  to  the  frequency  response 
parameters  such  as  bandwidth  or  amplitude  of  response  at 
resonant  frequency* 

In  the  analysis  of  a  servomeehanism.,  most  methods 
treat  the  servomeehanism  as  a  modification  of  a  system 
with  a  second  order  differential  equation.   The  output  of 
any  second  order  system  is  exactly  known  for  any  of  the 
useful  test  inputs  and  the  family  of  curves  is  not 
difficult  to  reproduce o   This  is  done  to  simplify  the 
description  of  the  output  with  a  given  input  to  one  that 
can  be  easily  formulated.   It  is  found  in  practice  that 
the  above  treatment  is  almost  always  a  useful  approximation 
if  it  is  but  remembered  that  it  is  an  approximation. 
5o2  The  Mathematical  Model  and  Graphs  of  Response 

An  early  task  of  the  designer  is  to  assemble  his 
likeliest  component ss    or  the  given  components <,  into  a 
representative  mathematical  model.   The  whole  process  of 
design  depends  upon  the  validity  of  the  model,,  being  only 
as  accurate  as  the  model  which  is  used  to  describe  the 
assembly  of  hardware.   The  transfer  functions  of  the 
components  can  be  gotten  from  experimental  tests  made  upon 
them  or  from  the  differential  equation  of  the  physics  of 
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the   component c     Even  when  the   form  of  the    differential 
equation  is  known,   experimentation  is   often  required   to 
obtain  the   values   of   the   parameters.      In  this  paper  <> 
analysis   is   described  as   the   determination  of   the    output/ 
input  characteristics   from  the   transfer  function,,   but  in  the 
determination  of   the   equations   for  a  component   all  of  the 
analysis   theory  and  methods  can  be   used  in  reverse o 

There  are  graphical  displays  which  are  used  to  portray 
either  the  mathematical  model  or  the  response  of  the  system 
to  a  given  input «>  Some  of  the  more  useful  ones  in  analysis 
and  synthesis  ares 

(a)  Frequency  Response   Graphs—    amplitude    of  the   output/* 
input  ratio  and  phase    shift  va.    frequency© 

(b)  Bode  Diagram  or  Attenuation  Diagram-«log  amplitude 
of  the  open  loop  output/input  ratio  and  phase  shift 
vs.    log  frequency o 

(c)  Nyquist  Diagram-  a  polar  plot  of  amplitude  of  the 
open  loop  output/input  ratio  and  phase  shift  with 
frequency  as   a  parameter c 

(d)  Inverse   Complex  Plane   Diagram- -a  polar  plot   of  the 
inverse    of  the    amplitude    of   the    open  loop   output/ 
input  ratio  and  phase    shift  with   frequency  as   a 
parametero 

(e)  Nichols  Chart-   a  log  amplitude   of  the   open  loop 
output/input  ratio  vs.   phase    angle   with  frequency 
as   a  parameter.      Log  amplitude    of  the   closed  loop 
(unity  feedback)   vs.   phase    shift   is   overlaid  on 
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the    chart o 
(f)      Root  Locus  Diagram--   the   complex  plane   upon  which 

the   poles   and  zeros   of  the   mathematical  model   are 

plottedo 
5c5      Specifications 

The  operation  of  a  servome onanism  can  be  described 
by  its  specifications.  These  specifications  are  defined 
in  such  a  way  as  to  describe  the  speed  of  response „  the 
stability,,  and  the  accuracy  of  the  servome onanism©  The 
specifications  used  in  Trial  and  Error  Design  Methods  all 
stem  from  the  ratio  of  the  output  of  the  servomechanism 
to   the    input  for  a  given  type   of  test  input© 

There    are    three   major  types   of   specifications e 
There   are    those    specifications  which  are   determined  from 
the   open-loop  frequency  responses    others  which  are 
determined   from  the   closed-loop   frequency  response  °s    and 
a  third   type   which  are   determined  from  the    transient 
response    to   a  specified  inputs   usually  a  unit   step 
function© 

Theoretically  the   closed-loop  response    to   a  step 

input  and  the   closed-loop   frequency  response   can  be    shown 

4 
to  be   equivalent      „      In  practice s    the   correlation  between 

the    two   is  usually  quite   remote g    and   the   conversion  from 

one    response    to  the    other  involves   graphical  methods   of 

integration  with  many  repeated  calculations o      There   is 

direct  correspondence   between  the    open-loop   and  the 
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closed=loop   frequency  response    and  this   relationship   is 
not  hard   to   determine,, 

There    are    a  great  many  different  measurements  used 
today  in  servome onanism  design   to   specify  the   character 
of  the    servo o      Most   of  these   measurements   are    inter- 
re  late  d5    and  many  of   them  are   either   synonymous   or  at 
least  they  describe   the    same   type   of  output* 

The    following  list   of   specifications   is   intended  to 
be   merely  a  sample    of  those  more   frequently  used* 
Many  of  them  require   amplifying  modifiers   not   to  be 
ambiguous. 

(a)    Open- Loop  Frequency  Response   Specifications 

(1)  Phase   Margin  is   180°  minus   the    angular 
difference   in  phase   between  the    output  and 
the    input  at  the   highest  frequency  where    the 
output  is   of  the    same    amplitude    as   the   input© 
It   is   a  measure    of  the   relative    stability  of 
the    system,,   being  unstable    at  non-positive 
values© 

(2)  Gain  Margin  is   the   ratio  of  the    amplitude 
of   the    input   to  the    amplitude    of  the    output 
at   that  frequency  where    the   phase    shift 
between  the    input  and   the    output   is   180° . 
It   is  usually  expressed   in  decibels   and  is 
a  measure    of  relative    stability,,    being 
unstable   for  non-positive  values   of  decibels© 


25 


odb 


9 

„-J?°      (1)  Phase  Margin 
(2)  Gain  Margin 


(b)  Closed-Loop  Frequency  Response  Specifications 

(3)  Peak  Amplitude  Ratio  (Mp)  is  the  maximum 
ratio  of  output  to  input 0   It  is  a  measure 
of  the  relative  stability  of  the  system* 

(4)  Frequency  of  Peak  Amplitude  Ratio  is  a 
measure  of  the  speed  of  response  of  the  systemo 

(5)  Bandwidth  is  one  of  the  specification  parameters 
which  has  no  standardized  definition<>   It  is  a 
measure  of  the  frequency  range  at  which  the 
amplitude  falls  within  specified  limits o   A 
popular  limit  is  that  the  amplitude  ratio  be 
between  4-  3db  and  =3db0   It  is  a  measure  of 
speed  of  response* 

(6)  Cutoff  Frequency  is  much  like  bandwidtho   It 
is  the  upper  frequency  at  which  the  amplitude 
ratio  reaches  some  specified  value o   Common 
ones  are  0  dbs  -3db5  ~6db5  or  -20db« 
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(3) >Tx  ^    Peak  ^P11*^6   Ratio 


-^^-  -  -f-   -\  (4)   Frequency  of  Peak 

.         \  Amplitude   Ratio 

(5)— J J  \ 

L      |    |   \     (5)  Bandwidth 

/^  ^ %r     (6)  Cutoff  Frequency 

( 7&8 )  Damping  Ratio  (  J )  and  Undamped  Resonant 
Frequency  (UJ%)    are  specification  factors 
in  that  they  completely  specify  the 
response  of  a  linear  second  order  servo- 
mechanism,.   All  of  the  other  specifications 
listed  are  fixed  by  (  f  )  and  (U)n)  for  a 
second  order  system.,   The  vast  majority  of 
servome onanisms  have  a  dominant  factor 
which  is  second  order  in  character* 
Consequently  ( f  )  and  {UJn)    specify  the 
type  of  response  of  a  higher  order  system 
to  a  more  or  less  close  degree,, 
(c)  Transient  Response  Specifications  to  a  Unit  Step 
Input 

(9)  Peak  Overshoot  (Mpt)  is  the  maximum  amplitude 
of  the  first  overshoot  measured  from  the 
final  steady  state  outputo   It  is  a  measure  of 
relative  stability °3    the  more  the  overshoot,  the 
less  stable  the  systemo   A  system  with  no  over- 
shoot is  said  to  be  "overdamped"  while  one  with 
an  overshoot  is  said  to  be  "underdampedM« 
(10)  Time  of  Peak  Overshoot  (t_.)  is  the  time  from 
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step  input  until  the  time  of  the  maximum 
amplitude  of  the  first  overshoot.   It  is  a 
measure  of  the  speed  of  response  of  the 
system. 

(11)  Rise  Time  is  the  time  for  the  output  of  the 
system  to  cancel  a  certain  portion  of  the 
error.   A  commonly  accepted  Rise  Time  is  the 
time  for  the  system  to  move  from  10$  of  the 
final  value  to  90$  of  the  final  value.  It 
is  a  measure  of  the  speed  of  response. 

(12)  Characteristic  time  ff  c  or  JfU)n)  is  the 
logarithmic  decrement  of  the  dominant 
portion  of  the  response.   It  is  a  measure 
of  speed  of  response. 

(13)  Settling  time  is  that  time  required  for  the 
servome onanism  to  reduce  the  error  below  and 
remain  less  than  a  certain  percentage  of  the 
input  step.  Commonly  accepted  percentages 
are  50    29   and  1  which  make  the  settling  time 
become  approximately  3S  45  and  5  times  the 
characteristic  time.   It  is  a  measure  of  the 
speed  of  response  of  the  system. 

(14)  Numbe r  of  Overshoots  is  the  number  of  over- 
shoots of  the  output  during  the  settling 
time.   It  is  a  measure  of  the  relative 
stability  of  the  system. 
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(15)    The  Error  Coefficient   is   a  measure   of  the 
steady  state   error  of  the    system.,      It   is 
dependent  upon  the   type    of  the    system  and 
the   gain  of  the   system., 
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3«4  Analysis 

A  very  important  part  of  the  synthesis  of  a  servo- 
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mechanism  is  the  analysis  of  the  system  under  test  to  see 
if  it  does  or  does  not  meet  the  required  specifications <, 
As  we  have  seen,,  these  specifications  are  invariably  an 
expression  of  the  comparison  of  one  form  of  output-to- 
input  relationship  or  another.   The  more  quickly  one  is 
able  to  determine  this  relationship^  the  shorter  the  job 
of  synthesis  will  be o   Many  of  the  analysis  methods 
available  are  in  reality  shorthand  techniques  at 
approximating  the  output- to- input  relationships o 
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Remember  that   linear   servome onanisms  have   output-to- 
input  relationships  which  are   expressible    as   linear 
differential  equations •      The    type    of  input   to  which  these 
systems   are   subjected  are   the   boundary  conditions  for  the 
differential  equations e      Since   it  has  been  found  that  the 
equations   can  be   handled  or  solved  more   readily  in  the 
Laplace   transformed  condition,   mainly  because   the   operators 
can  be  manipulated  algebraically,    the    system  equations  are 
normally  in  this   form*      This   is  known  as   the    frequency 
domain,,    and  it   is   the   one   in  which  most   synthesis   is  per- 
formed* 

The   frequency  response   of  a  system  is  easily  obtained 
by  simply  solving  the   closed-loop   transfer  function  for 
jCU*      If  the   open-loop   transfer  function  is  known  and 
available,,   rather  than  the   closed-loop   transfer  function* 
one    can  obtain  the   frequency  response    directly  without 
first  solving  for  the   closed-loop   transfer  function*      This 
is   done   by  plotting  the   frequency  response   of  the    open-loop 
system  and   then  transforming  the   coordinates*      If  there 
happens   to  be    a  feedback  function  other  than  unity,    the 
transformation  requires   one   more   step«    but  it   is   still 
useful* 

Unfortunately  many  of  the   specifications   are   pre- 
scribed in  the   time   domain  using  the    transient  response   to 
an  aperiodic   input.      This   situation  has   come    abqut   largely 
because    it   is  easier  for  the   specifications  writer,    and 
indeed  for  anyone,    to  visualize    the   effect  upon,    and  the 
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action  of,  the  servomechanism  in  that  plane*   So  there  is 
a  requirement  to  know  the  response  of  the  system  to  a 
transient  input.   There  are  several  avenues  of  approach 
available  to  determine  the  transient  response* 

The  most  obvious  approach 5  but  unfortunately  the  one 
which  usually  entails  the  most  labor,,  is  actually  to  solve 
the  differential  equation  with  the  prescribed  boundary 
conditions©   Normally  this  amounts  to  taking  the  Inverse 
Laplace  transform  of  the  closed-loop  transfer  function 
multiplied  by  the  Laplace  transform  of  the  input *  A 
modification  of  this  technique  is  to  get  the  approximate 
output  response  by  solving  only  for  the  dominant  features 
of  response 0    ignoring  the  less  important  features*  Note 
that  some  prior  experience  is  helpful  in  determining  what 
is  a  less  important  feature* 

A  second  avenue  of  approach  is  to  determine  first  the 
frequency  response  and  then  plot  the  time  response  from  the 
frequency  response*   It  has  been  shown  that  the  two 
responses  are  uniquely  equivalent*   It  is  extremely  un- 
fortunate that  the  equivalency  is  too  obscure  to  be  seen 
immediately  by  the  average  eye--or  for  that  matter  by  many 
a  trained  eye*   The  transformation  from  the  frequency 
response  to  the  transient  response  to  an  impulse  or  step 
involves  repeated  graphical  integrations 9    and  without  a 
digital  computer  can  be  a  tedious  operation* 

A  third p  and  the  most  frequently  used  approach  is  to 
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assume  that  the  system  is  like  a  second  order  system  and 
that  the  relationship  among  its  transfer  functions  its 
frequency  response^  and  its  transient  response  are  closely 
related  to  a  second  order  system.   The  great  majority  of 
servomechanisms  are  dominated  by  one  pair  of  complex 
conjugate  roots  possibly  with  a  single  dipole  near  the 
origin*   If  the  system  is  assumed  to  act  like  a  second  order 
one  with  only  the  complex  conjugate  roots ,  full  advantage 
can  be  taken  of  the  known  relationships  of  second  order 
systems  with  respect  to  transfer  functions „   frequency 
response  characteristics  and  transient  response  characterise 
tics.   The  results  can  be  appropriately  modified  if  there 
is  a  dipole  near  the  origin* 

The  frequency  response  can  be  determined  from  the 
transient  response*   This  is  of  value  when  a  physical 
component  is  available  for  transient  response  testings,  and 
the  transfer  function  is  unknown. 

In  summary.,  it  is  seen  that  the  loop  transfer  functions 
the  frequency  response^  and  the  transient  response  of  a 
servome chan ism  each  uniquely  defines  its  characteristics© 
Analysis  is  then  the  method  of  jumping  from  one  of  these 
four  descriptions  of  the  system  to  another. 
3  *  5  Synthe  sis 

The  determination  of  a  transfer  function  or  a  set  of 
components  which  will  comply  to  given  specifications  is 
synthesis*   Often  our  choice  of  components  is  either 
limited  or  partially  specified*   This  restricts  our  ability 
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to  determine  the  transfer  function  or  the  remaining 
components • 

In  synthesis  as  well  as  in  analysis9  heavy  reliance 
on  second  order  (or  possibly  3rd  order)  approximations  is 
required  in  many  of  the  present  methods,,  expecially  if  the 
specifications  are  given  in  terms  of  the  transient  response. 

An  analysis  of  the  specifications  will  quickly  show 
what  type  of  closed-loop  transfer  function  is  required  in 
order  to  remain  within  the  specifications-   The  real  trick 
of  synthesis  is  to  translate  this  closed-loop  transfer 
function*,  either  real  or  implied9  into  an  open-loop  function 
which  can  then  be  stated  in  terms  of  components* 

This  transformation  is  not  a  unique  one  and9  therefore,, 
there  is  an  infinity  of  combinations  of  components  which 
will  fall  within  the  specif icationso  Unfortunately  there 
is  a  much  greater  infinity  of  combinations  which  will  not 
fall  within  the  specif icationsc 

The  designer  is  allowed  a  considerable  amount  of 
leeway  in  the  selection  of  a  transfer  function  or  of  the 
components  he  must  usee   This  leeway  is  allowed  partly 
because  the  specifications  usually  require  the  system  to 
fall  within  certain  fairly  broad  limits  and  partly  because 
the  change  in  specifications  due  to  a  change  in  servo- 
mechanism  parameters  is  usually  not  a  sharp  change  <>   Since  P 
to  improve  upon  one  specif icationj  the  servome onanism  must 
often  be  allowed  to  relax  another  specification,,  and  since 
it  is  unlikely  that  the  specifications  can  be  described  as 
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optimum  to  begin  with,,  the  design  of  servome onanisms  is  as 
much  an  art  as  a  science «> 

The  most  usual  case  is  the  one  in  which  the  whole 
system,,  except  for  a  compensator  or  two,,  is  already  chosen,, 
and  it  is  only  a  matter  of  deciding  what  the  compensator  will 
bee  Even  the  selection  of  the  place  to  put  the  compensation 
may  be  restricted  due  to  the  limited  physical  accessability 
of  the  signal  flow. 
5o6  Compensation 

There  are  two  major  methods  of  compensating  a  servo» 
mechanismo   One  is  to  place  the  compensator  in  the  forward 
path  (called  cascade  compensation) 9    and  the  other  is  to 
place  the  compensator  in  the  feedback  path  (called  feedback 
compensation) 0 

Cascade  compensation  is  the  most  widely  used  of  the  two 
because  it  generally  requires  simpler  and  less  expensive 
elementso   It  is  easier  to  synthesize  a  system  using 
cascade  compensation  because  the  relationship  between  open= 
loop  and  closed~loop  frequency  responses  and  transfer 
functions  is  more  direct*  Drift  of  the  parameters  of  the 
active  elements  in  the  forward  path  disturbs  the  effect  of 
cascade  compensation  more  seriously  than  it  would  disturb 
feedback  compensation^©  Cascade  compensation  elements  are 
normally  placed  in  the  forward  path  at  signal  power  levels 
rather  than  at  output  power  levels e   This  allows  the  use  of 
smaller  network  elements  since  little  power  must  be 
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dissipatedo   The  two  types  of  cascade  compensation  that  are 
generally  used  are  lead  compensators  and  lag  compensators o 

Lead  compensators  act  as  a  derivative  path  in  parallel 
with  a  direct  patho   The  derivative  path  tends  to  increase 
the  output  to  input  ratio  at  high  frequencies..   At  these 
same  frequencies  it  decreases  the  amount  of  phase  shift 
that  would  occur  in  the  uncompensated  system©   The  impor= 
tance  of  the  lead  compensator  is  that  it  tends  to  cause  the 
180°  phase  shift  to  occur  at  a  higher  frequency  providing 
a  larger  degree  of  phase  margin  than  would  otherwise  be 
available  a   The  main  objection  to  lead  compensation  is 
that  it  causes  a  great  attenuation  of  the  forward  signal, 
requiring  a  substantial  increase  in  required  gain  in  the 
power  element  in  order  to  maintain  the  same  steady  state 
accuracy  as  before  compensation,.   Also  lead  compensators 
tend  to  increase  the  bandwidth  of  the  servome onanism 
making  it  more  susceptible  to  high  frequency  noise.   The 
transfer  function  of  the  compensator  is  of  the  forms 

Lag  compensators  act  as  an  integrating  path  in  par- 
allel with  a  direct  path*   The  integral  path  tends  to 
increase  the  output  to  input  ratio  at  low  frequencies 
while  the  phase  shift  is  increased  at  the  same  frequencies,. 
The  lag  compensator  is  used  to  increase  the  steady  state 
accuracy  of  the  uncompensated  system©   The  time  constant 
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of  the  polo  of  the  lag  compensator  must  be  made  large 
enough  that  the  accompanying  phase  shift  does  not  impair 
the  stability  of  the  system,,  and  yet  it  cannot  be  too 
large  or  the  transient  error  will  tail-off  too  slowly© 
The  transfer  function  of  the  compensator  is  of  the  form? 

Lead- lag  compensators  complement  the  virtues  of  the 
lead  compensator  with  those  of  the  lag  compensatoro   The 
relative  stability  is  increased  with  the  lead  compensator 
while  the  steady  state  accuracy  is  maintained,,  or  at  least 
the  attenuation  of  the  lead  compensator  is  offset 0   with  the 
lag  compensator* 

Combinations  of  lead  and  lag  compensators  can  be  used 
to  gain  further  benefits  from  these  compensators..   It  must 
be  remembered  that  isolation  must  exist  between  adjacent 
compensators  for  the  transfer  functions  to  be  correct  when 
multiplied©   This  isolation  can  be  achieved  with  a  buffer 
amplifier  or  cathode  followero 

Feedback  compensation^  while  more  complex  and 
expensive  in  components s   allows  a  greater  flexability  than 
cascade  compensation©   At  reasonably  large  gains,,  it  tends 
to  nullify  the  effect  of  the  forward  components  around  which 
it  is  placedp  making  the  whole  appear  like  the  inverse  of 
the  feedback  function..   Thus  it  is  an  ideal  method  of 
replacing  an  undesirable  function  with  a  desirable  one. 
A  forward  component  with  a  shifting  or  uncertain  gain  or 
parameters  in  the  transfer  function  can  be  made  quite 
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rigjid  by  the  proper  use  of  feedback  compensation,. 

Tachometer  feedback  compensation  is  a  very  common 
method  of  dampening  a  servomechanism  to  make  it  more  stable 
or  to  reduce  the  transient  overshoot*,   The  steady  state 
accuracy  suffers  with  tachometer  feedbacks  but  this  may  be 
corrected  by  placing  a  filter  in  cascade  with  the  tachom- 
eter in  the  feedback  loop* 
3  o  7  Summary 

The  Trial  and  Error  Design  Method  might  be  summarized 
as  follows? 

(a)  The  given  specifications  ares   input  signal 5 
desired  outputs  disturbances,,  allowable  error,,  plant 
elements ,  and  degree  of  stability  required© 

(b)  On  the  basis  of  experience  or  preliminary 
approximations,  select  a  form  of  compensation© 

(c)  Establish  the  parameters  for  the  compensation 
exclusive  of  system  gain© 

(d)  Adjust  the  system  gain  in  accordance  with  the 
stability  requirements o 

(e)  Analyze  system  to  see  if  the  error  is  satisfactory, 

(f)  If  the  error  exceeds  allowable  limitsj,  repeat  the 
process,,  using  different  compensation^  continue  until  error 
specifications  are  met© 

Although  the  theoretical  design  objective  is  the  attainment 
of  the  "best  possible"  servo,,  a  more  mundane  design  objec- 
tive is  the  attainment  of  an  "adequate"  servo© 
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CHAPTER  4 
BASIC  ANALYTICAL  DESIGN  METHODS 
4ol  General 

The  authors  received  almost  all  their  training  and 
schooling  in  servos  using  Trial  and  Error  Design  Methodso 
The  particular  problem  that  motivated  them  for  this  type 
of  thesis  was  the  fact  that  they  were  unable  to  recognize 
an  inconsistent  set  of  specificationse   Background  training 
was  more  than  ample  to  select  some  method  and  to  forge 
ahead  to  a  de signs  if  it  failed  to  satisfy  the  specifica- 
tion., to  start  overs  to  keep  starting  over  until  a  solu- 
tion was  found  or  patience  was  exhausted;,  never  really 
knowing  whether  a  solution  existed,,  or  in  the  case  where 
a  solution  was  found,,  whether  it  was  the  best  solution*. 

Analytical  Design  Methods  purport  to  solve  this 
dilemma o   And  this  is  true  in  the  sense  that  if  it  is 
assumed  that  a  performance  index  is  able  to  incorporate 
all  the  specifications  into  it5  the  method,,  being  a  pure 
mathematical  minimizing  or  maximizing  process,,  is  able  to 
immediately  reveal 5  solely  by  the  mathematics 0   whether  a 
stated  performance  index  can  be  achieved* 

Recall  that  the  three  main  themes  in  specifications 
are  speed  of  response 9    stability  and  accuracy!  it  seems 
reasonable  to  believe  that  some  performance  index  that 
includes  error D    time  and  an  appropriate  weighting  function 
for  the  error,  can  represent  the  specifications 0    or  at 
least  most  of  theme 
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In  selecting  a  performance  index 9   we  need  an  answer 
to  the  question  of  what  kind  of  output  is  desirable  for 
the  servoo   If  there  were  no  uncontrolled  disturbances 
(or  noise)  our  goal  would  be  to  make  the  output  follow  the 
input  perfectly,,   But,,  generally.,  in  the  presence  of 
disturbances;,  if  the  servo  follows  the  input  perfectly,,  it 
will  also  do  a  good  job  of  following  the  noise.   To 
establish  a  figure  of  merit  or  performance  index 9    then,, 
requires  a  compromise „   The  performance  index  must  be 
practicable  s   not  too  difficult  to  apply,,  and  of  general 
applicability,.  It  should  also  be  a  measure  of  the  average 
behavior  of  the  servo s    rather  than  be  affected  by  short~ 
lived  deviations  from  the  mean,,  or  shifts  in  the  time  axis, 
Quite  a  number  of  performance  indices  are  in  use  or  have 
been  proposed,,  some  of  which  are  discussed  below* 
4C 2  Performance  Indices  for  Transient  Signals 

(a)  Performance  Indices  which  are  not  time  weigh teds 

(1)  "Pi  -  f  ^(^)dtt  (called  Control  Area) 

(2)  p    [t^x  a/JSlX.  (called  Integral 

r2  -J0  Ie^Vdir       Absolute  Error  -  IAE) 

(3)  T^  _  C    (S>2/-fc)  J4-     (called  Integral- Square 

3    "J0  ^    wut  Error  -   ISE) 

Two   of  the    above    indices    (T-^  and  P3)    can  be  used   in  purely 
analytical  procedures  whereas   P2  contains  discontinuities • 
?2_  was  proposed  by  T„   Mo    Stout75    Pg  by  Pickeison  &  Stout8,, 
and  P-   by  Hall     and  Sartor iuso      P3   is   sometimes   called 
the  Hall-Sartorius   criterion©      It  can  be    seen  that   all  of 

these   indices  heavily  weight   the   error  at  the   beginning  of 
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the  transient.  Even  though  they  favor  rapid  speed  of 
response,,  they  also  tend  to  cause  adjustments  that  are9  in 
general^,  less  damped  than  is  usually  requiredo   From  this 
viewpoint  an  IAB  index  is  better  than  an  ISE  indexo   IAE 
is3  however,,  normally  treated  with  the  use  of  analog 
compute rsc  Note  that  none  of  these  indices B    though,,  place 
a  penalty  on  errors  occurring  late  in  the  transient,,  and 
thus 3  there  is  no  good  assurance  that  the  servo  does  not 
have  a  long  tail  offo   To  provide  for  thisp  and  at  the  same 
time j,  recognizing  that  no  physical  system  can  respond 
instantaneously,,  another  set  of  performance  indices  are 
time  weightedo 

(b)  Performance  Indices  which  are  time  weighted^ 


(4)  Bt  ■=  J0"t  e&)d±.  (called  Time  Weighted 

Control  Area) 

(5)  "PL-  J  -fc  I^WI  J-fc     (called  Integral- Time- 

°  Multiplied  Absolute 

*<*>    ^  Error-  ITAE) 

(6)  K  ~   J0"^  ^  tjd't         (called  Time  Weighted 

ISE) 

11  1  o 

P4  was  proposed  by  Nims   and  P&  by  Graham  and  Lathrop   , 
Againp  P4  and  Pg  can  be  treated  in  purely  analytical 
procedures;,  whereas  Pg  is  usually  treated  with  an  analog 
computer*.   Pg  can  also  be  treated  with  a  set  of  what  are 
called  standard  forms  of  which  the  foremost  proponents 
are  Graham  and  Lathrop*   It  Is  discussed  further  in 
Chapter  6o   The  major  objection  to  these  indices  is  that 
they  tend  to  tip  the  scales  the  other  ways   they  put  too 
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much  weight  on  the  latter  portion  of  the  transient©   They 
might  be  more  practical,,  if  they  ignored  error  after  it 
reached  a  certain  limito   Unfortunately  such  an  expression 
would  again  not  be  analytic alo   However,,  P5  is  quite 
popular  in  the  literature  for  computer  use. 

(c)  Another  index  proposed  by  Nims1   iss 

However D   we  can  find  no  other  literature g   for  or  against  itc 
It  can  be  used  where  the  objective  is  to  make  P^  and  P4 
some  value  other  than  zero.  P-  (ISE)  is  by  far  the  most 
popular  index  for  purely  analytical  methods.   One  reason 
for  this  is  that  error  becomes  more  undesirable  as  It 
increases  in  magnitude •  As  the  below  figure  shows,, 
weighting  the  square  of  the  error  rapidly  penalizes  large 


error. 


4 


Square-  crtrcr 


HOD-  s<?U4r<L~ 

ELrrov* 

However 0   probably  the   greatest  reason  for  the  wide  usage 

of  ISE    (and  Mean  Square  Error  discussed  in  the   next 

paragraph)    is  mathematical  convenience.      There    is   a  highly 

polished  body  of  mathematical  knowledge   that  has  been 

developed  particularly  around   the    idea  of  a  mean  square 

value  o 
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4c5  Performance  Indices  for  Stochastic  Signals 

It  was  indicated  earlier  that  usual  specifications 
concern  the  response  of  the  system  to  typical  inputs p 
whereas  systems  are  actually  subjected  to  random  inputs. 
Thi's  is  one  of  the  basic  reasons  for  interest  in  statisti- 
cal criteria,.   By  far  the  most  commonly  used  and  described 

index  is  the  root~me  an- square  criterion,,  first  proposed 

9 
by  A,  Co  Hall  and  developed  by  James 0   Nichols  and 

Phillips5 o   It  is  characterized  by  its  mean  square  value „ 
thuslys  r  T 

(8) 


Hidd-b 


Another  criterion  suggested  by  Oldenburg  and  Sartorius15 
is  characterized  by  the  following  integral^ 

The  use  of  the  root-mean- square  criterion  was  in- 
spired by  the  writings  of  Wiener14  and  it  has  been  highly 
developed  by  many  writers?  of  whom  two  good  ones  are 
Newton2  and  Laning15* 

Recall  that  determination  of  the  response  of  a  system 

to  a  stochastic  input  does  not  yield  a  function  of  time5 

but  rather  a  statistical  characterization  of  the  output 

signal  ensemble o   Stochastic  signals  consist  of  two 

classes? 

(1)  The  process  is  stationary  if  the  statistical 

behavior  of  the   process   that  generates   the   ensemble    is 

independent   of   time* 
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(2)  The  process  is  non~ stationary  if  the  statistical 
behavior  varies  with  timeo 

Many  books  are  written  on  the  fact  that  if  a  stochastic 
process  is  stationary^  it  is  possible  to  use  a  single  member 
of  the  ensemble  of  the  process  to  determine  the  process 
statistics.   (Various  books  also  justify  the  assumption  that 
stationary  signals  are  valid  for  use  in  servo  work.   No  one 
has  successfully  "cracked  the  barrier"  of  non- stationary 
signals  yets  except  for  special  cases.)  Furthermore s    the 
average  value  of  a  signal  from  a  stationary  process  can  be 
found  either  by  taking  an  ensemble  average  at  a  particular 
time  or  by  taking  the  time  average  of  a  single  member  of  the 
ensemble.   Because  of  this  (see  Appendix  B) g    the  mean-square 
value  of  a  stationary  signal  is  precisely  equal  to  its  auto- 
correlation function,,  evaluated  with  the  argument  equal  to 
zero.   If  Pg  is  then  selected  as  the  performance  index 9    an 
incredibly  powerful  tool  results. 

The  rms  value  of  a  random  quantity  can  often  be 
calculated  practically  under  conditions  in  which  many  of  its 
other  statistical  properties  cannot.  For  example  <,  if  a 
random  quantity  has  a  normal  probability  distribution,,  all 
statistical  properties  are  determined  from  its  rms  and  mean 
value . 

It  should  be  stressed,,  thought  that  the  rms  value  of  the 
error  does  not  characterize  the  error  complete ly\,  and  use  of 
this  criteria  may  result  in  overlooking  some  important 
aspects  of  the  problem,   e  (t)  is  independent  of  the 
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distribution  of  e(t)  in  the  frequency  spectrum.   Since  a 
system  usually  operates  with  other  systems  connected  to  its 
the  other  systems  will  transmit  e(t)  in  a  manner  that  depends 
on  their  dominating  frequency  and  on  the  spectrum  of  e(t). 
Gille-1-6  indicates  as  an  example  the  stabilization  of  an 
airplane  carrying  passengers*   Resonant  frequencies  of  the 
order  of  magnitude  of  the  duration  of  a  human  pace  should  be 
avoided  because  such  frequencies  are  the  most  adverse  to 
comforto   It  is  just  as  important  to  consider  the  frequency 
spectrum  of  the  error  as  its  rms  value.   Nevertheless^,  there 
is  widespread  use  of  PQ  in  design  analysis.   This  seems  to 
be  due  less  to  its  intrinsic  worth  as  a  criterion  than  to  the 
convenience  attached  to  its  calculation. 
4.3  Configurations 

By  establishing  standard  configurations-,  much  of  the 
mathematics  can  be  worked  out  for  one  time  and  only  the 
results  are  then  necessary  for  design.   The  configurations 
used  ares 

(a)  Fixed  Configurations 

In  this  configuration;,  all  the  physical  elements , 
including  the  compensation,,  are  essentially  fixed.   Op- 
timization procedure  then  consists  of  the  adjustment  of  a 
few  free  parameters  (such  as  gain,,  or  one  time  constant  of  a 
filter,,  or  perhaps  the  ratio  of  a  pair  of  time  constants  of 
a  filter)  so  as  to  minimize  the  performance  index.   This 
technique  is  frequently  used  in  conjunction  with  the  Trial 
and  Error  Methods  where  the  form  of  the  system  has  been 
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fixed  by  other  considerations  and  only  the  best  numerical 
values  of  the  free  parameters  are  sought o   On  the  surface 
it  appears  easy  to  use  since  the  only  requirement  is  to 
differentiate  the  performance  index  with  respect  to  the 
free^  parameters  and  to  set  the  partial  derivatives  equal 
to  zero.   This  technique  is  also  easy  to  apply  in  connec- 
tion with  the  analog  computer*) 

(b)  Free  Configuration 

In  this  configuration,,  quite  a  high  degree  of 
mathematical  sophistication  is  called  for,,  because  here 
the  entire  transfer  function  (or  weighting  function,,  if  in 
the  time  domain)  of  the  system  is  allowed  to  vary  in 
minimizing  the  performance  index •   This  can  only  be  done  by 
the  use  of  the  calculus  of  variations «   This  is  one  of 
Wiener's  big  contributions ,  producing  an  implicit  transfer 
function  in  the  form  of  an  integral  equation  known  as  the 
Wiener-Hopf  equation,.   Solution  of  this  equation  is 
extremely  difficulty  Wiener  conceived  a  process  for  solution 
called  "spectrum  factorization".   Although  this  kind  of 
procedure  is  very  satisfying  intellectually;,  it9  in 
general,,  leads  to  a  transfer  function  which  is  physically 
unrealizable «, 

(c)  Semi-Free   Configuration 

This  can  be  thought  of  as  the  usual  type  design 
problem,,  and  is  the  type  that  one  first  thinks  of  in  the 
Trial  and  Error  Design  vernacular;   All  the  plant  elements 
are  specified  as  fixed  elements,,  but  there  is  no  constraint 
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on  the  type  of  compensation  allowed.   This  technique  again 
frequently  leads  to  a  compensation  which  is  physically 
unrealizable p  but  it  has  a  very  real  practical  advantages 
it  provides  a  goal  to  aim  for  in  compensating  by  Trial  and 
Error  Methods*   In  providing  the  goal,,  it  frequently 
indicates  the  direction  one  must  take  in  the  selection  of  a 
compensator  type© 
4«>4   Rudiments  of  Analytical  Design 

All  of  the  methods  use  the  overall  or  closed-loop 
transfer  function.,  Adopt  the  following  symbology  (taken 
from  Newton  )  (Same  as  that  on  symbols  page.,  but  repeated 
here  for  convenience)  g    v^j  /<A  .  (Jj(v) 
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Because  the  closed  loop  function  is  used,  it  sometimes  is 
not  too  easy  to  manipulate  back  to  the  original  block 
diagram.   However „    this  disadvantage  (if,  in  factp  it  is  one) 
is  offset  by  the  fact  that  these  methods  do  not  involve 
solving  for  the  roots  of  the  characteristic  equation,  unless 
they  are  desired  in  the  final  completed  paper  design.   This 
is  a  very  real  advantage  if  the  characteristic  equation  gets 
above  an  order  of  about  five. 

Practically  speaking,,  the  application  of  these  methods 
can  lead  to  failure  if  the  limitations  of  the  mathematical 
model  are  not  continuously  kept  in  mindo   They  almost 
inevitably  lead  to  non-linear  operation  because  the  per- 
formance  index  demands  minimum  error o   The  minimization  of 
dynamic  error  calls  for  higher  gain,  leading  to  saturations 
and  it  calls  for  cancellation,  or  near  cancellation  of  the 
plant  elements  by  the  compensator  elements B    leading  to  wide 
bandwidtho   (For  example  <,  design  of  a  second  order  system  by 
these  techniques  will  invariably  lead  to  infinite  gain  which 
is  not  only  physically  impossible,  but  practically.,  undesir- 
able o  ) 

(a)  Transient  Input 

Let  us   first  consider  procedures   for  a  transient 
input*     Use   of  the   ISE  criterion  has  been  highly  developed 
by  Newton  &  Gould2*4 2   use   of  ITAE  criterion  has  been 
developed  by  Graham  &  Lathrop12    (discussed  in  Chapter  6)0 
The   ISE  criterion  is  represented  ass 
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By  Parseval^s   Theorem,,    the    Integral  can  be   rewritten  ass 
T      =    J-      C*~e(s)1e(-S)ds 

where  Ys(s)    is   the  Fourier  transform  of  y9  (t).      If  Ya(s)    is 
a  rational   function,   Iy  can  be   written  in  the   forms 

where  C(s)    and  D(s)    are   polynomials   in  s0     Definite   Integrals 
of  this   form  have   been  evaluated  and  tabulated  in  terms   of 
the   coefficients   of  the   polynomials    (for  example s    see  Appendix 
E   of  reference   2)»     Evaluation  of  the    integral  then  givess 

x3  «  ^  ft,  H  - . .  .yu\ 

where   P  represents   the   free  parameter© 

Iy   can  be   minimized  by  solving  the    set  of  K  equations   of  the 

forms  —\  -p 

Except  for  the  simplest  forms 9    solution  of  this  set  of 
equations  can  be  difficult,,   Many  times-,  the  easiest  approach 
is  simply  to  plot  Iy  versus  Pk5  holding  other  parameters 
constants,  the  minimum  value  can  be  determined  closely  enough 
for  engineering  purpose s«   An  important  condition  for  use  of 
the  developed  integral  tables  is  that  D(s)  must  have  all  its 
zeros  (roots)  in  the  left  half  plane •   This  is  also  the 
condition  for  a  stable  system*   Therefore 9    the  parameters 
must  not  be  allowed  to  take  on  such  values  as  to  cause  D(s) 
to  have  zeros  In  the  right  hand  plane*   Otherwise  the  integral 
table  Is  invalid  and,,  besides;,  the  system  is  unstable*, 
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Note  that  another  limitation  exists  above?   Ye(s)  must 
be  a  rational  function*   (It  is  not  meant  to  cite  this  as  a 
disadvantage  compared  to  Trial  and  Error  Methods^  since 
those  methods j,  in  general.,  require  rational  functions  also©) 
See  Chapter  69   under  Newton7 s  Method,,  for  a  way  around  this 
limitation 

The  above  scheme  can  only  be  used  with  fixed  configura~ 
tion  systems o  For  semi-free  configurations  with  transient 
inputs 0  Newton  uses  translation  functions »   This  is  also 
postponed  to  Chapter  6« 

(b)  Stochastic  Input 

Let  v(t)  and  q(t)  represent  stochastic  input  and 
output  respectively  of  a  linear  system  whose  weighting 
function  is  w(t)«   Prom  the  convolution  integral  we  know  that 

From  the  definition  of  autocorrelation  functions  (Appendix  B) , 
and  from  manipulating  the  above  integrals 9    the  following 
relation  is  obtained? 

which  expresses  the  autocorrelation  function  of  the  output 
in  terms  of  the  autocorrelation  function  of  the  input  and  the 
weighting  function* 

In  a  similar  manners 


t<A   =  I 
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#**)  =J^^)^v(^-t,)d^ 


which  gives  cross- correlation  between  the    input   and  the   out- 
put in  terms   of  autocorrelation  of  the   input   and  the 
weighting  function,.      If  the   correlation  functions   and  the 
weighting  functions   are   Fourier  transformable,    then 

The    latter  expression  is   frequently  used   to  evaluate  W(s), 
since    it   is   the   ratio  of   the    eross-power-density  spectrum  to 
input-power-density  spectrumo 

Prom  the   block  diagram,  tje  (t)    ~   i(t)«^.(t)o      The  mean- 
square-error  is   identically  the  value   of  the    autocorrelation 
function  of  error,   with'T  3   Os 

After   suitable  manipulation,    the   power  density  spectrum  is 
obtained  ass 

This  equation  is    suitable   for  many  variations «      For  example, 
if  'Lr&)-  ^SW+^nt^ where  "Md($  represents   the    data  or   signal 
component   and  UJ^(-Q    the   noise   component,    and  we   assume   them 
to  be  uncorrelated,    thens 

(Obtained  by  letting  §>,„$  =^j(i)+^W$|and  §J$=$ii($*&ui(%) 

in  the   previous  equations,, ) 

Now,    the   problem,    similar  to   that  described  for   transient 

input  is   to   find  the    area  beneath  the   error-power- density 
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spectrum  along  the  imaginary  axis  and  then  minimize  this 
error.   This  can  be  done  by  evaluating  and  minimizing  the 
following  integral; 


Of  course j,  the  complexity  of  the  mathematics  involved  has 
increased  greatly  over  that  for  the  fixed  configuration 
transient  signal „   Although  the  error  may  again  be  minimized 
by  setting  the  partial  derivatives  equal  to  zero5  it  is 
usually  best  to  proceed  with  a  series  of  plots  of  error 
versus  the  free  parameters,. 

(c)  Stochastic  Input o  Free  and  Semi-free  Configurations 
Newton  shows  that  the  actual  output  can  be  eliminated 
from  the  error  expression,  and  the  error  then  is  a  function 
only  of  the  input  and  the  desired  output  as  follows 8 

He  then  proceeds  to  prove  that  minimization  of  the  error  for 
a  physically  realizable  weighting  function  requires  sat° 
isf action  of  this  expressions 

This   is  the  Wiener-Hopf  integral  equation  in  the    time   domain*, 
of  which  UJ^(t)    is   the    implicit   solution5   where  Uj^(t) 
represents   that  weighting  function  that  minimizes    U^  (^tj  » 

It  is   an  understatement  to  say  that   solution  of  this 
equation  is   difficult c      There    are   a  few  rare    situations 
where    its   solution  in  the    time   domain  is   self«evidento 
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Newton  (Chapter  5,  reference  2)  has  a  very  good  expose  on 
the  explicit  solution  to  this  equation  in  the  frequency 
domain.   Suffice  to  say  that  its  solution  consists  of 
Wiener3 s  "spectrum  factorization"  technique  whereby  a 
function  is  split  in  such  a  manner  that  one  portion  consists 
of  that  component  whose  poles  lie  in  the  LHP  and  the  other 
portion,  that  which  lies  in  RHP.   By  discarding  the  RHP 
portions,  the  explicit  solution  of  Y\Jm(s)  gives  a 
guaranteed  stable  transfer  function,,   This  solution  in  the 
frequency  domain  is  much  easier  if  the  functions  are  Fourier 
transformable*.   Some  authors  have  developed  a  solution 
whereby  the  whole  derivation  for  optimum  transfer  function  is 
conducted  in  the  frequency  domain. 

Since  solution  for  a  free  configuration  is  somewhat 
academic,  the  equation  has  been  modified  to  provide  for  the 
semi-free  configuration  We  quote  that  solution  he  res 

f     &(-s)£v;fej    7 

Although  this  expression  is  quite  ominous  at  first  glance, 
a  breakdown  of  the  symbols  helpss 

Wcm(s)  -  optimum  cascade  compensation  for  minimizing  U^Uz) 

Csr      "  fixed  or  plant  elements 
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/\   (s)    ~   any  factor  of  Z\.(s)  which  includes   all   the   poles 

and  zeros   of  /\    (s)  in  the   LHP. 

O(^)^."   component  of    O    (s)  which  has   all   its  poles   in 

LHP   such  that  Q    (s)  -     O    (s)^  has   all  its  poles 
in  the   RHP. 

ZVfe)~    A(s>)/4    ($J  s   remaining  factor  of  A(s) 

which  includes  poles   and  zeros   of  £J>    (s)    in  RHP0 
So s    if  the   correlation  functions   are   known ^    then  solution  of 
the   above   equation  for  Wcm(s)    becomes   a  mathematical 
factoring  problem,,   where    stability  is   guaranteed. 

If  the   fixed  elements   are  minimum  phase   functions   the 
above   expression  reduces   tog 


<Lfr\ 


(s)    . 


where  Wm(s)   =  1  ^yC  fe)/<gv*(s)  J+  /<&/•  (V 

For  our  purposes  a  minimum  phase  transfer  function  is  one 
which  has  no  zeros  in  the  RHPo   Now5  Wm(s)  is  the  solution 


for  the  overall  system  transfer  function  for  minimum  U^C^J 
where  there  are  no  fixed  elements e   Therefore 9    the  overall 
system  transfer  function  is  independent  of  G^(s)  if  G^(s) 
is  minimum  phase.  What  this  means  is  that  when  the  desired 
ouptut  of  a  system  is  equal  to  the  inputs  the  above  equation 
will  always  call  for  a  Wcm(s)  such  that  the  overall  system 
transfer  function  is  1  (thus  giving  a  me an- square -error  of 
zero)  if  the  fixed  elements  are  minimum  phase.   And  this 
sounds  reasonable j,  since  the  ultimate  goal  of  any  compensa- 
tion is  to  cancel  all  attenuation  and  phase  shift  caused  by 
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the  plant  over  the  whole  frequency  spectrum,  in  order  that 
the  output  will  be  an  exact  replica  of  the  input.  Even 
though  this  is  the  ultimate  goal,,  wes  of  course  recognize 
that  physically^  the  goal  is  unattainable.   Hence 0    this 
method  in  general  leads  to  a  compensating  function  which  is 
physically  unattainable.   The  above  methods  rather  non- 
chalantly assumed  that  the  correlation  functions  were 
available  j,  whereas  actually.*  the  obtaining  of  a  correlation 
function  in  analytical  form  is  another  barrier  to  be 
surmounte  d , 
4,4  Correlation  Functions 

(a)  Derivation  from  Theoretical  Considerationsg 
Appendix  B  lists  some  autocorrelation  functions  that 
can  be  derived.   Because  of  this$,  they  are  popular 
for  use  in  design*  Furthermore 3  their  use  does 

not  represent  too  great  a  departure  from  the  world 
of  practicality, 

(b)  Derivation  from  Experiments 

Any  sophisticated  or  complicated  design  calls  for 
the  determination  of  the  correlation  function  by 
experiment,   Newton**  presents  an  approximate 
numerical  procedure  for  computing  them  from 
oscillograph  traces.   He  also  indicates  an  analog 
computer  method  whereby  a  stylus  is  manually 
moved  over  the  oscillograph  tracings.   He  states 
that  currently  the  M, I,T,  Servome onanisms 
Laboratory  prefers  a  numerical  procedure  because 
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computers  can  bo  use do 

After  determination  by  experiment „  the 
functions  are  frequently  not  analytical  and  some 
approximation  scheme  is  necessary  if  they  are  to 
be  used  in  analytical  equations 0   For  example 9 
reference  17  describes  a  scheme  for  approximating 
a  correlation  function  by  a  series  of  damped 
cosine  functions •   If  a  function  is  approximated 
by  any  such  scheme 9   generally  the  whole  design 
problem  is  of  such  size  that  computers  are 
necessary  for  the  entire  design  problem© 

4  c  5   Summary 

The  Analytical  Design  Method  might  be  summarized  as 

follows s 

(a)  The  given  specifications  ares   input  signal,, 
desired  output,,  disturbances,  performance  index 
and  required  value  of  same,  plant  elements,,  and 
degree  of  freedom  allowed  in  compensationo 

(b)  Classify  the  problem  according  to  free,  semi- 
free  s  or  fixed  configuration.. 

(c)  For  free  or  semi-free  configuration,,  use  an 
appropriately  derived  formula,,  For  fixed 
configuration,  express  the  performance  index 
as  a  function  of  the  free  parameters!  minimize 
the  performance  index  by  adjusting  the  param- 
eters o 
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(d)  See  if  the  compensation  thus  determined  yields  the 
required  value  of  performance  index •   If  it  does 
not.,  the  specifications  cannot  be  met?  if  it  does, 
practical  realization  may  begin© 
The  design  objective  is  the  attainment  of  the  "best 
possible"  servo  to  meet  the  specifications,  using  the  chosen 
performance  index  as  the  criterion. 
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CHAPTER  5 
TRIAL  AND  ERROR  DESIGN  METHOD  REFINEMENTS 

5d  Introduction 

The  purpose  of  Chapter  5  Is  to  aid  a  serv ©mechanisms 
designer  to  find  an  appropriate  method  of  approaching  the 
solution  of  his  particular  design  problem.   To  this  end,,  a 
summarization  of  the  techniques  found  profitable  by  the 
authors  is  made  along  with  references  to  more  complete  study 
of  the  technique  involvedo 

The  design  problem,,  when  determined  by  the  de signer g 
will  contain  three  partss   (1)  the  specifications  which 
must  be  metj,  (2)  fixed  elements  to  be  included,,  and  (3) 
any  limitations  or  preferences  in  the  type  of  compensators 
to  be  used*  The  open-loop  transfer  function  or  frequency 
response  of  the  simplest  possible  servo  containing  the 
required  elements  should  be  analyzed  immediately  in  order 
to  decide  the  amount  and  type  of  compensation  required© 
This  is  done  by  synthesizing  a  system  which  will  meet  the 
specifications  and  then,,  using  the  compensators  and 
variable  gains 9    constrain  the  actual  servo  model  to  be 
similar  to  the  system  first  synthesized,.   The  final,,  or 
possibly  feedback^  step  is  to  analyze  the  final  servo 
model  to  insure  that  it  does  indeed  meet  the  specifica- 
tions o   More  care  must  be  taken  with  approximations  in  this 
final  analysis  than  in  the  preliminary  analyse 3. 

Following  a  guide  to  design  methods  available 9    some 
of  the  less  well  publicized  methods  are  discussedo   Some 
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root  locus  tochniquosp  Mitrovic^s  method,,  and  some  trans- 
formation  methods  are  explained  to  the  extent  that  they  can 
be  used  with  the  help  of  this  paper©   Theory  and  develop- 
ment of  the  methods  is  left  to  the  references  lis  ted  j,  with 
only  a  hint  as  to  the  validity  of  the  methods<>   A 
knowledge  of  the  basic  frequency  domain  and  time  domain 
methods  is  assumed*  The  best  general  review  of  the  methods 
available  is  either  "Ordnance  Engineering  Design  Handbook4" 
or  "Handbook  of  Automation  Computation  and  Control 9 
Volume  I19" 
5«2  Guide  to  Design 

The  various  methods  to  design  a  servomechanism  can  be 
listed  in  any  number  of  ways?  type  transfer  function  used,, 
domain  used,  the  graphical  plane  used.,  analysis  or 
synthesis 3  etc©   None  of  these  groupings  is  an  aid  to  the 
designer  when  he  is  attempting  to  make  his  choice  for  a 
design  problem  at  hand*.   The  breakdown  of  methods.,  shown 
on  the  following  pagesp  while  redundant,  does  allow  the 
designer  to  concentrate  his  attention  to  the  methods  found 
fruitful  by  the  authors©   It  is  by  no  means  an  exhaustive 
collection* 
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AB1 
AB2 

AB3 
AE1 

AE2 

AF1 
AP2 


Me  thod 

Reference 

root  locus 

16  ,,25 

solve  equation 
G 

1+5 

3 

Mitrovic 

pare  5«5 

Bode  diagram 

3 

Nyquist 
criterion 

20 

Bode  diagram 

3 

solve  equation 
GKjw) 

20 

Limitations 


unity  feedback 
unity  feedback 


BA1 

root  locus 

pare  5«4 

unity  feedback 

BA2 

solve  equation 
1-W 

3 

unity  feedback 

BC1 

root  locus 

par*.  5.3 

BC2 

second  order 

approximation 

3 

BC3 

third  order 

approxlmati  on 

par<>  5*3 
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Mo  thod 

Reference 

Limitations 

BC4 

inverse  Laplace 
transform 

20 

BD1 

Bode  diagram 

3 

BD2 

root  locus 

par*  5«3 

BD3 

second  order 
approximation 

3,16 

BD4 

solve  equation 
W(Jw) 

20 

BE1 

root  locus 

16 

BE2 

Routh  criteria 

19  9    20 

CB1 

Bode  diagram 

3 

CB2 

second  order 
approximation 

3,  4,  16,  1 

CB3 

curve  fitting 

20 

CD1 

Nixon's  method 

par«  5«7(a) 

CD2 

Chestnut  &  Mayer 

pare  5o7(b) 

DB1 

Bode  diagram 

3 

DB2 

second  order 
approximation 

3,  19 

DG1 

Floyd0 s  method 

par. 

5, 

►  6(a) 

DC  2 

Guillemin5  s 

par. 

5, 

»6(b) 

DC3 

Stallard 

4 

DF1 

Nicholas  chart 

3 

must  know 
feedback 

DF2 

Nyquist  diagram 

20 

must  know 
feedback 

DP3 

Inverse  polar  plane 
plot 

20 

must  know 
feedback 
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Me  thod 

Reference 

Limitations 

PA1 

Bode  diagram 

3 

PCI 

Chestnut  &  Mayer 
diagrams 

1*  4,  19 

FD1 

Nicholas  chart 

3 

FD2 

Nyquist  diagram 

20 

FD3 

Inverse  polar  plan 

9 

plot  20 

(Notes   paragraphs  indicated  above  refer  to  this  paper.) 

5»5  Root  Locus  Analysis 

The  root  locus  method  was  first  developed  by  Walter 
Evans-1-0  in  1948„   It  has  gained  considerable  popularity 
with  usage.   There  have  been  many  modifications  of,,  and 
techniques  developed  about 9    the  root  locus • 

Basic  uses  of  the  root  locus  are  the  plotting  of  the 
open-loop  transfer  function  poles  and  zeros  on  the  complex 
plane  and  the  plotting  of  the  locus  of  the  roots  of  the 
closed-loop  transfer  function  on  the  complex  plane  using 
the  poles  and  zeros  of  the  open-loop  transfer  function* 

The  plotting  of  the  root  locus  is  basic  to  the  general 

3  19  20 
method  and  is  assumed  to  be  known  by  the  designer  '  9        • 

The  limits  of  gain  which  will  allow  a  system  to 
remain  stable 9    and  the  second  order  approximation  of  the 
relative  stability  for  a  given  gain,  can  be  determined 
from  the  root  locus. 

Analysis  of  the  root  locus.,  from  the  position  of  the 
closed-loop  poles  and  zeros,  can  be  determined  by  the 
following  methods  which  are  listed  approximately  in  order 
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of  increasing  difficulty,  but  also  in  order  of  increasing 
accuracy* 

(a)  Graphical  Representation  of  Dominant  Pole-Zero 
Locations  Versus  Selected  Specificationsc 

Elgerd  and  Stephens21  have  graphed  the  step  input 
response  for  almost  every  combination  of  roots  and  zeros  up 
to  fourth  order  for  selected  distances  from  the  origin  and 
for  selected  conjugate  pole  distances© 

Abbott  and  Pat ton22  have  graphed  the  time  for  the  out- 
put to  first  equal  the  inputs  the  peak  overshoot,,  the  time 
of  peak  overshoot ,  and  the  settling  time  (all  for  a  step 
input)  for  closed-loop  systems  containing  a  complex 
conjugate  root  pair  (1)  alone  (second  order)  (2)  with  a 
zero  (second  order),  or  (3)  a  real  root  (third  order) o 
They  further  show  that  systems  containing  a  complex 
conjugate  root  pair,  a  real  root,  and  one  or  two  zeros  can 
be  broken  into  partial  fractions  of  the  above  three  types 
and  the  transient  responses  then  added. 

Both  of  the  above  techniques  require  the  system  to 
be  normalized  to  make  one  of  the  roots  or  zeros  unity 0 
Elgerd  and  Stephens  consider  more  types  of  systems  and 
their  article  (Trans  AIEE,  Vol  78,  Pt  II,  1959)  is  more 
apt  to  be  available,  at  this  time.   The  method  of  Abbott 
and  Patton  requires  essentially  no  interpolation  due  to 
their  mapping  of  the  responses  on  the  complex  plane e   Both 
methods  show,  along  with  Chu's  (see  below)  that  a  real 
pole  or  zero  loses  most  of  its  dominant  influence  when 
it  is  more  than  twice  as  far  from  the  origin  than  one  of 
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the  complex  conjugate  roots •   See  Figure  1  at  the 

end  of  this  chapter. 

23 
(b)  Chu's  Approximation  Equations 

This  method  computes,  either  graphically  or  analyt- 
ically j,  the  dominant  terms  of  the  transient  response  to  a 
step  input o   Prom  this  he  finds  simple  equations  for  the 
peak  overshoot  and  the  time  of  peak  overshoot.  Let 
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Aio) 


z   a  factorable  constant 

=  the  product  of  the  distances  from  all 

the  zeros  to  the  origin 
~  the  product  of  the  distance  from  all  the 

roots  to  the  origin 

"  G"o+.  tf^Oo  the  location  of  one  of  the 
dominant  complex  roots 
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A„ 


S0)   s  the  product  of  the  distances  from  all 

of  the  zeros  to  S6 
i 

^CSa)       ~   the  product  of  the  distances  from  all 
the  other  roots  to  Sa 

V<So)      =  A*jA(So)  -Ancj  B^soj  -Ang  Sa 

^"C&o)  ~    ^**    sum  °^  *ke    angles   from  the   real 

axis   to  the   line   from  the   zeros   to  s    . 
ritlQ  OfeJ-   the    sum  of  the    angles   from  the   real 

axis   to  the    line    from  the   other  roots 


to   s0, 


f\fi4   <£*   -  the  angle  from  the  positive  real  axis 
to  the  line  from  the  origin  to  s0„ 
/1(^\   ^  the  product  of  the  distances  from  all 
the  zeros  to  J7 
D  ^-j  =  the  product  of  the  distances  from  all 
the  other  roots  to  G^   0 
(c)  Wheeler" s  Approximation  Equations24 

Wheeler's  equations  are  similar  to  Chu?Sj,  except 
that  he  makes  the  substitutions 

K  "     Acs.)  Ac^ 

|<  .AisL    s  C^^  ■  Fi**l  value- 

Thus  j,    if  one   knows   the   final  value    of     CftJ  or  F^C $'&>   he 

may  plot   the    step  response    from  the   ope n=  loop  poles   and 

the   closed- loop   roots o      In  any  casep   he    can  plot  the   error 

to   a  stepc 
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whores   K^»ft\  -  the  product  of  the  distances  for  all  the  open- 
loop  poles  to  Sb« 
(d)  Graphical  Residue25 

The  graphical  residue  gives  the  true  transient 
response  to  an  input  step*   If  there  are  repeated  rootss  it 
will  not  work  in  the  form  given  below  but  the  above  reference 
should  be  consulted! 
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(e)  Frequency  Response  from  Root  Locus 

The  frequency  response  of  the  system  can  be  gotten 
graphically  for  any  given  Ui  by  measuring  the  distances  from 
all  the  roots  and  zeros  to  jw  on  the  imaginary  axisc 
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5*4   Root  Locus  Synthesis 

Having  found  where  the  desirable  locations  for  the 
dominant  roots  are ,  the  designer  must  now  be  able  to 
"maneuver"  the  system's  roots  into  these  favorable  locations* 
This  can  be  attempted  by  varying  the  parameters  of  the  given 
system  or  by  adding  a  compensator  in  cascade « 

(a)  Effect  of  Varying  One  Parameter  of  Root  Locus24 
Dr.  Wheeler  demonstrated  that  any  one  parameter  of 

the  open-loop  equation  may  be  factored  out  of  part  of  the 
characteristic  equation  and  then  be  used  as  the  gain  of  the 
final  root  locus t    thus  showing  the  effect  on  the  system  of 
varying  this  parameter*  After  factoring  out  the  variable 
parameter,  a  root  locus  plotter  would  be  a  definite  aid* 
This  type  of  plotter  is  available  in  different  degrees  of 
sophistication26*   The  real  use  of  this  technique  would  be  to 
see  the  true  effect  of  one  variable  parameter* 

(b)  Placement  of  Lead  or  Lag  Compensators2"7 

The  Ross-Warren  technique  is  a  natural  outgrowth  of 
the  root  locus  equations  and  of  earlier  efforts  by  Walters28 
and  Aseltine29*   The  locus  is  found  to  a  point  where  the 
dominant  complex  conjugate  roots  are  desired  and  the  phase 
angle  is  noted.   The  placement  of  the  lead  or  lag  elements 
is  partially  fixed  by  this  phase  angle  because  they  must  be 
placed  so  as  to  change  that  angle  to  -180°.   The  compensator 
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is   further.)    and  completely,    fixed  by  stipulating  that   the 
error  coefficient    (Kv)    be  unchanged  by  the   compensation* 
This   stipulation  is  unnecessarily  restrictive   and  is   one   of 
the   major  drawbacks   of  the   methodo 

Measure:       £)Cj     (jCj  {jj^  )    Cosf'/CJ    K, 

if   -  Un-i)TT-£)c  ujhev^.  n   is   such  +ha+  |<PI  <  lQo° 


2    =    OJn     K  sin/} /6C  ^/'n© 


whe  re  s 

C?c        =   ^e   phase   angle   measured   to   the   desired 
root   location0 

CzIq^       s   the   product   of  distances   from  the   zeros 
to  the   desired  root  location  divided  by 
the   product   of   the   distances   from  the 
poles   to  the   desired  root   location., 

^ci  Jc  ~   &e3Cr>1-I>tlon  of  tne    desired  root,    in 

second  order   terminology* 
rx  ~   the   gain  of  the   uncompensated  equation* 

Z?  ~   zero   location  of  the   compensator 

-   pole    location  of  the   compensator 

3?  /p  °   attenuation  of  the   compensator* 

A  modification  of  the    above   method   allows   the 

designer  to   approximate    immediately,    a  lead  filter  with  a 
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zero j,  and  to  use  other  than  the  uncompensated  Ky*   See 
Figure  2  at  the  end  of  this  chapter  with  which  one  enters 
with  the  angle  to  be  compensated  and  obtains  the  zero 
placement^  .   The  effect  of  a  real  pole  can  be  approximated 
from  the  same  graph  by  subtracting  the  angle  found  at  the 
pole  distance  from  the  angle  found  at  the  zero  distance. 

Mariotti  has  developed  a  similar  method™  for  the 
determination  of  pole -zero  location*   The  Ross-Warren 
technique 5  modified,,  requires  successive  approximations  to 
place  the  compensator  at  the  place  with  the  best  error 
coef ficient,,  while  Mariotti  has  graphed  this  placement  so 
that  it  can  be  accomplished  immediately«>   His  work  shows 
that  the  error  coefficient  increases  as  the  compensator  is 
moved  further  from  the  imaginary  axis. 
5o5  The  Method  of  Mitrovlc5*31 

(a)  Theory 

The  method  of  synthesis  presented  by  Dusan 
Mitrovic  is  a  different  and  potentially  powerful  approach  to 
the  solution  of  the  synthesis  problem*   This  method  utilizes 
the  characteristic  equation  in  the  unfactored  polynomial  forme 

Separating  the  real  and  the  imaginary  parts  of  the 
characteristic  equation  and  solving  for  a0  and  a-^  produces 
two  parametric  equations  in  terms  of  the  higher  order 
coefficients  and  with  the  variables  j  and  (jQ  e   Inspection 
of  these  equations  reveals  that  each  power  of  (A)  contains  a 
single  coefficient  of  the  characteristic  equation  and  an 

68 


easily  tabulated  function  of  T   • 

The  basic  consideration  of  the  Mitrovic  method  is  to 
let  a0  and  a^  become  variable  and  to  observe  the  (a^sa0) 
plane  0 

(b)  Properties 

The  locus  f  ( ifj  of  the  parametric  equations  is 
plotted  for  a  specific  Jf  and  is  found  to  have  some 
remarkable  properties s 

(i)   r(o) 

The  servome onanism  is  stable  if  certain  conditions 
are  met.  This  is  because  the  total  phase  rotation 
as  CO  -*  00  is  fixed  for  a  stable  system  similar  in 
nature  to  the  Nyquist  criteria*  This  curve 
immediately  shows  the  range  of  a-j_  and  a0  for  which 
the  system  will  be  stable,, 

<2>   r(fk) 

The   characteristic  equation  has  no  roots  with  J 
greater  than   fy  if  certain  conditions   are  met. 
This   curve    immediately  shows   the   range   of  a-^   and 
aQ  for  which  the    system  will  have   a  J    less   than 
J/^     •      More    important j,    it  will   show   the   exact 
a^   and  a0  necessary  to  make    the    system  have   roots 
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at   any  given  (JJ  n  and      J   ^» 

(3)  ro) 

The  characteristic  equation  has  all  real  roots  if 
certain  conditions  are  me to   Here  again 9    the  range 
of  a^  and  aQ  are  shown.   The  real  roots  are 
quickly  determinable.   The  effect  of  changing 
al  0I*  ao  3-s  iMnediately  evident*   If  there  is 
but  one  complex  pole  pair9  it  can  be  determined 
easily  once  the  real  roots  are  known. 
(4)   Third  order  equations 

If  the  system  can  be  approximated  by  a  third 
order  equation  (i.e.  a  dominant  complex  conjugate 
pole  pair  and  a  maximum  of  one  dominant  real  root 
with  no  more  than  3  dominant  zeros)  this  method 
is  considered  the  ultimate  in  flexibility  by  the 
authors » 

A  third  order  characteristic  equation  can  be 
normalized  so  that  a^  and  &2   are  both  one 0   Now 
all  third  order  equations  can  be  solved  by  one 
set  of  I  (Jy°   This  set  is  tabulated  once  and 
for  all3531o   The  ax  and  a0  required  for  all 
dominant  root  combinations  is  immediately 
discernable.   Thus  a  very  powerful  tool  is 
available  to  determine  the  parameters  of  any 
system  with  a  closed-loop  characteristic  equation 
of  third  order  or  less. 
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(c)  Plotting  of  PUT) 

The  equations  developed  in  the  theory  can  be  used  to 
plot  CCS)  with  a  table  of  $*($)* 

An  easier  method  is  available e    thougha  using  a 
table  that  has  the  ^ftCf)  UJ    already  constructed 
for  the  Jf     in  question.   Tables  II „   III9  IV  and  V 
(at  the  end  of  this  chapter)  are  constructed  for 
values  of  Jf  of  0,  0o5P  0.7  and  loO. 

a,  =  a^.(Xw)  +  a3^co)-t--«+^nr^^ 

These  tables  are  constructed  for  values  of  LU 
from  zero  to  one.   To  insure  that  the  area  of 
interest  falls  within  the  above  iX)5    the  characteris- 
tic equation  must  be  normalized  so  that  both  an 
and  an=^  are  unity.   This  is  accomplished  easily 
by  the  substitution  s  s  an=i  PT9  and  it  insures 
(for  equations  with  all  negative  real  parts  in  the 
roots)  that  all  roots  lie  between  0  and  =lo   The 
transformed  parameter  is  related  to  the  actual 
parameter  bys 
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(d)    Determining  Stability 

The  fT(d\ curve   is  plotted  and  point    (a^pa0)    is 
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noted. 

The    system  is   stable   ifs 

(1)  The   point    (al9a0)    lies   in  the   first 
quadrant  9   and 

(2)  while (JJ  varies  between  Q  and  ^00  the 
curve   alternately  cuts   the    lines   aQ  and 

a^  with  the   provision  that  a     is  cut  firsts 
and 

(3)  that  the  total  number  of  points  of  inter- 
section is  nP  where  n  is  the  order  of  the 
characteristic  equation. 

(e)    Determining   Jj 

The  r~(jfrjcurve   is  plotted  and  the   point    (a^9a0)    is 
noted.      All  roots  have   damping  factor  greater  than    Jj^  ifs 

(1)  The   point    (aljr    a0)    lies   in  the    first 
quadrant „    and 

(2)  while   curve  r(jicjalternately  cuts   the    lines 
aQ  and  a^  with  the   provision  that  a0   is  cut 
firsts    and 

(3)  that   the    total  number  of  points   of  intersec- 
tion is  m5   where  m  is   the    largest   integer 
which  fulfills   the    inequality? 


m  + 


[,_(-irj§-<n(t+  4^) 


where   n  is   the   order  of  the   characteristic 

equation  and      Q  i   -af*    —  ^OS    JK  • 

If  the   point    (a11Ja0)    lies   on  thefu/ycurvep 
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the  damping  factor  of  a  complex  conjugate 
root  pair  is  Jfy     and  LU^   is  the  fjj   of  that 
point  on  the  curve  0 

(f)  Determining  the  Real  Roots 

The  nCy  curve  is  plotted  and  the  point  (a^5aQ)  is 
note do 

All  real  roots  lie  on  the  |  (')  curve  at  a  distance  0**2 

(when  JP  °  i.o,  Ul  z  <T~  )    from  the  origin  such 
that  a  line  tangent  to  &X   will  pass  through  (al3a0). 
Conversely 5  at  any  place  that  a  line  from  (a^9a0) 
is  tangent  to  the  curve 0  there  is  a  root  on  the 
negative  real  axis  at  v      0 

(g)  Lead  and  Lag  Compensation 

The  use  of  the  Mitrovic  method  is  quite  easy  with 
lead  or  lag  compensators •   The  characteristic  equation 
of  a  lag  compensated  transfer  function  will  probably 
have  the  lag  pole  in  coefficients  other  than  a^  or  a0„ 
It  will  be  found  that  the  effect  of  the  lag  pole  on 
these  coefficients  is  so  small  that  it  can  be  deletedo 

The  characteristic  equation  of  a  lead  compensated 
system  has  the  same  problem  of  having  the  compensator 
pole  in  the  higher  order  coefficients.)   Unfortunate ly9 
it  is  a  larger  term,  and  cannot  be  ignored*   Instead  it 
must  be  fixed,,  in  order  to  get  numerical  values  for 
the  coefficients o   This  would  also  have  to  be  done  to 
the  zeros  and  the  gain  if  they  appeared  in  the  higher 
order  terms.   Then  any  computations  involving 
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variations  of  parameters  which  were  fixed  in  the 
higher  order  terms  would  either  be  approximations  or 
the  curve  would  have  to  be  continuously  replotted* 

5o6  Frequency  Response  to  Transient  Response 

The  several  methods  for  obtaining  the  transient 

response  to  an  impulse  input  from  the  frequency  response 

curve  make  use  of  the  following  relationships 

f  (i)  =    M«*  ujf  &  [.<*<&«*»)]  i"J 

whe re  ^ni.LF(^w) J      is   the   real  part  of  the   frequency  response© 

The   real  part   of  the    frequency  response    can  be   computed 
directly  from  the   closed=loop   frequency  response    curve*      If 
the   frequency  response   were   plotted  in  polar  coordinates,, 
the   real  part  could  be   read  directly*      Chen  and  Shen^ 
suggest   a  chart  which  has   the    closed-loop   frequency 
response   plotted  in  polar  coordinates   and   overlaid  upon  it 
is   the   open-loop p   unity  feedback  frequency  response© 
G-ould*   recommends   a  charts    similar  to   a  Nichols   chart,,    on 
which  the    open-loop   frequency  response   curve    is  plotted  as 
log-gain  vso   phase    angle    (See   Figure   3   at  the   end  of  this 
chapter.)      Overlaid  upon  the   chart   is   the   magnitude    of  the 
real  part   of  the    closed-loop   frequency  response    for  a 

unity  feedback  system* 

4   19 
(a)    Floyd? s  Method 

The   real  part  of  the   closed-loop  frequency 

response   curve   is   approximated  by  the    sum  of  a  series   of 

trapezoids? 
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f(*)  =  If  ^ 


€*/ 


sin  LQ^fe 


gig 


4i± 

4l* 


2. 


A  check  on  the  correctness  of  the  trapezoids  iss 

k 


-f(o)  =  -f  ^Trua^ 


(b)    Guillemin's  Method  4*25 

The    straight   line    approximation  of  the   real  part 
of  the   closed-loop  frequency  response   is   differentiated 
twice  s    leaving   a  series   of  impulses   of  height   CLl  and  at 
frequency    OJ £   s 


m  =  - 


7Ti' 


7. 


4-  ya^cosujjt 


A  check  on  the  correctness  of  the  impulses  is; 


(c)  A  Computer  Method^ 

Levadi  has   combined  Floyd's  method  and  Guillemin8 s 
method  in  a  way  that   is   easy  to  program  on  a  compute r« 
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The  output  from  the  computer  is  the  impulse  response,  the 
step  response ,  and  the  response  to  a  ramp  input • 
5»7  Transient  Response  to  Frequency  Response 
(a)  Nixon *s  Method20 

A  straight  line  approximation  of  the  transient 
response  to  a  step  input  is  made  and  it  is  broken  into  the 
sum  of  a  series  of  ramp  functions*   These  functions  could 
be  expressed  in  the  time  domain  as  a  magnitude  multiplied 
by  a  unit  ramp  with  a  time  delay*   Instead,,  they  are 
expressed  directly  in  the  frequency  domain  (by  Laplace 
transform)  and  multiplied  by  s.   This  is  the  transfer 
function  because s   before  multiplying  by  39    this  sum  is 
the  Laplace  equation  for  the  output  to  a  unit  step  input. 
If  the  transient  response  to  an  impulse  were  usedj,  the 
multiplication  by  s  would  have  been  unnecessary*   To  find 
the  frequency  response s    substitute  j w   for  ss 


ra«o  =  TL[&i£ 


Since    the  equation  is    indeterminate   at  CO  =   0P    LJHospital*s 
rule   must  be  used, 

(b)    Chestnut  and  Mayer vs  Method4 

The    transient  response   to   a  step   input   is  broken 
up   into   the    sum  of  a  series   of   step   functions  each  occurring 
at   an  equal   interval   of  times   but   being  of  a  variable 
height    ACl°o 
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1-=.! 


The  transient  response  to  an  Impulse  is  broken  up 
into  the  sum  of  a  series  of  pulses  of  uniform  width   At 
but  of  variable  height  C^s 
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TABLE   II 


TABLE  FOR  COMPUTING  MITROVIC'S    P£q)    CURVE 


U) 

w2 

# 

ft 

% 

% 

a 

t 

1ti 

% 

%  ' 

fio 

0 

0 

1 

0 

0 

0 

.05 

.003 

-1 

0 

.003 

0 

.10 

.01 

-1 

0 

.010 

0 

0 

.15 

.023 

-1 

0 

.023 

0 

7001 

.20 

.04 

-1 

0 

.040 

0 

T002 

0 

.25 

.063 

-1 

0 

.063 

0 

:004 

0 

0 

.30 

.09 

=  1 

0 

.090 

0 

7008 

0 

.001 

0 

.35 

.123 

-1 

0 

.123 

0 

7015 

0 

.002 

0 

0 

.40 

.16 

-1 

0 

.160 

0 

-026 

0 

.004 

0 

7001 

0 

.45 

.203 

-I 

0 

.203 

0 

7041 

0 

.008 

0 

7002 

0 

.50 

.25 

-1 

0 

•  250 

0 

7063 

0 

.016 

0 

7004 

0 

.55 

.303 

-1 

0 

.303 

0 

7092 

0 

.028 

0 

7008 

0 

.60 

.36 

-1 

0 

.360 

0 

7130 

0 

.047 

0 

7017 

0 

.65 

.423 

-1 

0 

.423 

0 

7179 

0 

.075 

0 

7032 

0 

.70 

.49 

-1 

0 

.490 

0 

7240 

0 

.118 

0 

7058 

0 

.75 

.563 

-1 

0 

.563 

0 

7316 

0 

.178 

0 

7100 

0 

.80 

.64 

-1 

0 

.640 

0 

7410 

0 

.262 

0 

7168 

0 

.85 

.723 

-1 

0 

.723 

0 

7522 

0 

.377 

0 

7273 

0 

.90 

.81 

-1 

0 

.810 

0 

7  656 

0 

.531 

0 

7431 

0 

.95 

.903 

-1 

0 

.903 

0 

"815 

0 

.735 

0 

7663 

0 

1.0 

1.0 

-1 

0 

1.0 

o- 

1.0 

0 

1.0 

0 

-1.0 

0 

U)    VS.    I//     forf^O 


*%(%*>)  =  *Mflu) 


k-i 


0,  -G*ii.OfJu)+Qi'H(J>u>)+  ~.  +  Q,tih($u)) 
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TABLE  III 


TABLE   FOR  COMPUTING  MITROVIC ■ S 


r(o.s) 


CURVE 


U) 

w2 

•Vx 

^2 

% 

14 

Y5 

^6 

% 

^8 

f* 

fio 

0 

0 

-1 

0 

.05 

.003 

„1 

.050 

0 

0 

•  10 

.01 

-1 

.100 

0 

7001 

0 

.15 

.023 

=  1 

.150 

0 

1003 

.001 

•  20 

.04 

=  1 

.200 

0 

1008 

.002 

.25 

.063 

-1 

.250 

0 

7016 

.004 

0 

0 

•  30 

.09 

-1 

.300 

0 

7027 

.008 

0 

7001 

.35 

.123 

-1 

.350 

0 

1043 

.015 

0 

1002 

.001 

.40 

.16 

-1 

.400 

0 

-064 

.026 

0 

7004 

.002 

0 

0 

.45 

.203 

-1 

.450 

0 

1091 

.041 

0 

1008 

.004 

0 

7001 

.50 

.25 

-1 

.500 

0 

1125 

.063 

0 

7016 

.008 

0 

7002 

.55 

.303 

-1 

.550 

0 

1166 

.092 

0 

■5028 

.015 

0 

1005 

.60 

.36 

-1 

.600 

0 

1216 

.130 

0 

7047 

.028 

0 

7010 

.65 

.423 

-1 

.650 

0 

7275 

.179 

0 

7Q75 

.049 

0 

1021 

.70 

.49 

-1 

.700 

0 

•2343 

.240 

0 

7118 

.082 

0 

1040 

.75 

.563 

-1 

.750 

0 

7422 

.316 

0 

7178 

.134 

0 

7075 

•  80 

.64 

-1 

.800 

0 

7512 

.410 

0 

7262 

.210 

0 

7134 

•  85 

.723 

-1 

•  850 

0 

•  614 

.522 

0 

7377 

.321 

0 

7232 

•  90 

.81 

-1 

•  900 

0 

7729 

•  656 

0 

7531 

•  478 

0 

7387 

•  95 

.903 

-1 

.950 

0 

■5  857 

.815 

0 

7735 

.698 

0 

7  630 

1.0 

1.0 

-1 

le00 

0 

1.0 

1.0 

0 

1.0 

i 1 

1.0 

0 

-1,0 

ti>    VS.  if   -for   j^O.S 

a,  =  a  act^+Q*  a  (t  w>  •  •  •  + Q*  •#•  (£  *>) 
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TABLE   IV 


TABLE  FOR   COMPUTING  MITROVIC'S 


r(c.i) 


CURVE 


u) 

w2 

fi 

^2 

V* 

t* 

Ys 

Y« 

fr 

a 

f. 

0.0 

.000 

4 

0 

0 

.05 

.003 

-1 

.070 

:048 

.10 

.010 

-1 

.140 

T096 

0 

.15 

.023 

-1 

.210 

1144 

.001 

.20 

.040 

-1 

.280 

1192 

0 

.002 

0 

.25 

.063 

-1 

.350 

^240 

7001 

.004 

7001 

0 

.30 

.090 

-1 

.420 

1288 

7002 

.008 

7003 

.001 

.35 

.123 

-1 

.490 

1336 

7002 

.016 

1007 

.002 

.40 

.160 

-1 

.560 

1384 

1004 

.027 

1014 

.004 

0 

.45 

.203 

-1 

.630 

1432 

1005 

.042 

?026 

.008 

0 

1002 

.50 

.250 

-1 

.700 

7480 

1007 

.065 

7044 

.014 

.001 

1004 

.55 

.303 

=  1 

.770 

7528 

7010 

.095 

7070 

.025 

.002 

1009 

.60 

.360 

-1 

.840 

-576 

7012 

.134 

7109 

.043 

.003 

1018 

.65 

.423 

-1 

.910 

7624 

7015 

.185 

7162 

.069 

.005 

■g034 

.70 

.490 

-1 

.980 

7  672 

1019 

.249 

7235 

.108 

.009 

1062 

.75 

.563 

-1 

1.050 

1720 

1024 

.328 

7332 

.163 

.015 

7108 

.80 

.640 

-1 

1.120 

7768 

1029 

.424 

7458 

.241 

.023 

7180 

.85 

.723 

-1 

1.190 

1816 

7034 

.541 

1^20 

.346 

.036 

o  CiPQ 

.90 

.810 

-1 

1.260 

7864 

7041 

.680 

825 

.488 

.054 

1463 

.95 

.903 

-1 

1.330 

7912 

7048 

.844 

-1.082 

.675 

.078 

7713 

1.0 

1.0 

-1 

1.400 

7960 

1056 

1.036 

~lc598 

.918 

.112 

i  - 

-L075 

0)    V*.  T//      far  f=^7 


-^0»  =  <pk(/)cu 


k-i 
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TABLE  V 


TABLE   FOR   COMPUTING  MITROVIC'S     r(l.o)     CURVE 


60 

fl 

y2 

n 

* 

ft 

"% 

W 

% 

~W 

0 

-1 

0 

0 

.05 

=1 

.100 

"008 

.000 

1000 

.10 

-1 

.200 

1030 

.004 

-.001 

.000 

.15 

-1 

.300 

^068 

.014 

1003 

.000 

.20 

-1 

.400 

1120 

.032 

1001 

.000 

.000 

.25 

=1 

.500 

-.188 

.062 

1020 

.006 

-sOOl 

.000 

.30 

=1 

.600 

-.270 

.108 

1041 

.015 

-S005 

.002 

.000 

.35 

=1 

.700 

".368 

.172 

1075 

.032 

?013 

.005 

-.002 

.40 

-1 

.800 

-.480 

.256 

1128 

.061 

<s028 

.013 

-.006 

.45 

-1 

.900 

-.608 

.364 

1205 

.111 

<s058 

.030 

-.013 

.50 

-1 

1.000 

1750 

.500 

1312 

.188 

•S109 

.062 

-.035 

.55 

=1 

1.100 

1908 

.666 

1458 

.302 

•S194 

.122 

-.076 

.60 

-1 

1.200 

-1.080 

.864 

1648 

.467 

<?327 

.224 

-.151 

.65 

1.300 

-1.268 

1.108 

1893 

.696 

"sOSo 

.392 

-.287 

.70 

=-1 

1.400 

-1.470 

1.372 

-1.201 

1.009 

*©  BZw 

.659 

-.518 

.75 

1.500 

-1.688 

1.688 

-1.582 

1.424 

-1.246 

1.068 

-.901 

.80 

-1 

1.600 

=1.92 

2.048 

-2.048 

1.966 

=1.835 

1.678 

-1.510 

.85 

1.700 

-2.168 

2.456 

-2.610 

2.662 

=2.640 

2.565  -2.453 

.90 

-1 

1.800 

-2.430 

2.916 

-3.281 

3.540 

=3.720 

3.826 

*3  oo  «  & 

.  y  O 

1.900 

-2.708 

3.430 

-4.073 

4.643 

-5.146 

5.586 

-5.971 

1.0 

-1 

2.000 

»3.0 

4.0 

-5.0 

6.0 

=  7.0 

8.0 

-9.0 

a.  = 


(JO     VS.  If     -for  j=I.O 


Q.  *  Q^.(^cuj+a5fs(^a))^...+a^(j;ai) 
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CHAPTER  6 
ANALYTICAL  DESIGN  METHOD  REFINEMENTS 
60 1  Introduction 

In  this  chapter  we  summarize  selected  methods  from 
current  literature  that  concern  Analytical  Design.   Some  of 
the  methods  do  not  strictly  meet  the  definitions  of  the 
Analytical  Design  Method  at  the  beginning  of  this  paper 9 
but  they  are  included  heres  because  the  literature 
generally  considers  any  technique  which  uses  a  performance 
index  as  its  main  criterion  to  be  "analytical". 
6.2  Newton7 s  Methods2 

Chapter  4  was  built  almost  exclusively  around  the 
writings  of  Newton,  Gould  and  Kaiser2*   However 9    they  have 
done  a  great  deal  more  than  to  develop  some  basic 
mathematics?  they  have  expended  much  effort  toward  applying 
the  basic  mathematics.   Their  methods  use  ISE  and  MSE 
exclusively,  since „  as  they  state „   these  indices  are  the 
only  ones  of  engineering  usefulness  that  lead  to  reasonably 
straightforward  mathematical  analysis. 
(a)  Translation  Functions 

Newton  prove s9   by  means  of  variational  calculus „ 
that  translation  functions  can  be  treated  in  a  manner 
similar  to  correlation  functions.   Translation  functions 
(see  appendix  B)  are  used  with  transient  inputs.   The 
formulas  indicated  in  Chapter  4,  then,  can  be  used  by 
merely  substituting  the  correct  translation  function  in 
the  place  of  the  correlation  function.   Most  translation 
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functions  can  be  found  by  Parseval's  Theorem©   Use  of  step 
and  ramp  inputs 9   even  though  they  have  no  Fourier  transforms 
impose  no  special  limitation,,  since  it  is  possible  to 
introduce  a  convergence  factor  to  find  a  Fourier  transform* 

(The  Fourier  transform  found  this  "way  is  the  same  as  the 
Laplace  transform;  however.,  special  care  is  required  in 
taking  the  inverse  transform).   This  scheme  is  also  much 
easier  to  use  if  the  fixed  elements  are  minimum  phase. 
What  appears  to  be  a  particular  advantage  of  the  method  is 
its  indifference  to  the  kind  of  transfer  function  that  exists 
in  the  plant*   For  example ,  it  is  able  to  handle  a 

Gjf  (s)sr  £"°  S  as  readily  as  a  «£f(s>)  -  ST  . .   Although  this 
method  normally  employs  the  frequency  domain  for  finding 
the  best  compensations  it  frequently  finds  the  minimum 
value  of  ISE  by  working  exclusively  in  the  time  domain. 
(b)  Saturation 

As  indicated  earlier*  compensation  by  error 
minimization  will  almost  invariably  lead  to  saturation  in 
a  linear  systems  rendering  the  assumed  mathematical  model 
invalid.   The  reason  for  this  is  that  the  nc ailed  for" 
compensation  transfer  function  is  usually  that  which  will 
cancel  the  effect  of  the  plant  transfer  function.   Since 
the  plant  normally  consists  of  integrations  and  lags,  this 
means  that  the  compensation  (or,  at  least 9    the  equivalent 
cascade  compensation)  will  call  for  several  differentiations. 
Differentiations  lead  to  large  signal  excursions  in  the 
systems  with  the  result  that  saturation  occurs  somewhere. 
Differentiation  also  calls  for  more  bandwidth,  since  the 
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signals  have  high  frequency  components o   Many  physical 
systems  have  the  characteristic  of  linearity  for  large 
amplitudes  at  low  frequencies,  but  can  only  tolerate  small 
amplitudes  at  higher  frequencies  if  linearity  is  to  be 
preserved.   To  make  the  Analytical  Design  Method  practical, 
this  reference  introduces  constraints  into  the  equations  to 
insure  that  saturation  does  not  occur. 

(l)  Transient  Inputs  -  The  most  direct  approach  for 
avoiding  saturation  is  to  limit  the  peak  value  of 
that  signal  in  the  mathematical  model  that 
corresponds  to  the  signal  in  the  physical  system 
that  is  likely  to  cause  saturation.   However,  it 
is  impractical  to  express  the  peak  value  of  a 
signal  as  a  function  of  the  free  parameters  for 
systems  above  second  order*   The  reason  is  that 
finding  the  roots  of  the  characteristic  equation 
above  that  order  requires  the  use  of  numerical 
values  so  specific  values  of  the  parameters  are 
necessary*   But  the  integral- square  value  of  a 
signal  can  frequently  be  expressed  in  terms  of  the 
free  parameters!  since  the  integral- square  value 
gives  large  weight  to  large  values^  then  some 
degree  of  control  of  the  peak  value  is  obtained* 
One  procedure  is  to  express  the  integral- square 
error  as  a  function  of  the  free  parameters °9    also 
express  the  integral- square  value  of  the  saturaration 
signal  as  a  function  of  the  free  parameters.   Then 
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solve  the  equations  by  adjustment  of  the  parameters. 
For  fixed  configurations  this  is  best  done  by  trial 
and  error.   For  free  and  semi-free  configurations « 
this  reference  introduces  a  more  refined  technique., 
employing  a  method  developed  by  LaGrange.   By  use  of 
a  synthetic  function  and  a  constant  called  the 
LaGrange  multiplier,  the  minimization  of  error 
problem  is  converted  from  one  with  a  constraint  to 
one  without  a  constraint .  Either  of  these  methods 
works  well  in  a  fixed  configuration  problem,,   such 
as  a  positional  servo.   In  this  type  of  servo,, 
saturation  most  frequently  results  from  demanding 
an  acceleration  of  the  output  member  or  load,  which 
the  motor  is  simply  not  capable  of  providing.   It 
is  to  be  noted  that  each  introduction  of  a  constraint 
effectively  eliminates  one  free  parameter  in  system 
design. 
(2)  Stochastic  Inputs  -  Newton" s  method  for  handling 
saturation  with  stochastic  inputs  is  quite 
impressive.   The  method  involves  the  derivation  of 
a  modified  Wiener-Hopf  equation  with  the  insertion 
of  a  LaGrange  multiplier.   Since  explicit  solution 
of  the  multiplier  is  difficult 9    graphical  procedures 
are  usually  resorted  to9  for  determination  of  the 
multiplier.  The  method  seems  to  have  no  restrictions 
as  long  as  the  use  of  the  rms  value  of  the 
saturating  signal  is  an  acceptable  criterion  for 
constraining  saturation.   If  this  is  unacceptable 9 
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then  there  seems  to  be  no  choice  <,    analytically 
speaking j,  but  to  separate  the  design  problem  into 
two  parts?  a  linear  and  a  non- linear  mathematical 
model,  employing  non-linear  techniques  after 
saturation  occurs • 
(c)  Bandwidth  Minimization 

The  above  paragraph  indicated  one  reason  why  design  by 
this  method  called  for  increased  bandwidth*   There  is 
another  reason?   Since  the  input  signals  are  represented  by 
their  power-density  spectra.,  or  energy- density  spectra 
(for  transient  inputs) p   minimization  of  ISE  or  MSE  will 
demand  that  the  system  transmit  all  those  frequencies  in 
the  input  with  negligible  amplitude  and  phase  change,   If 
the  input  signal  is  aperiodic ,  its  spectrum  includes  all 
frequencies 5  if  it  is  periodic,  it  includes  discrete 
frequencies  which  are  integer  multiples  of  the  fundamental 
out  towards  infinity,,   Hence  ,  for  minimum  error,  namely 
zero,  an  infinite  bandwidth  is  called  for* 

Excessive  bandwidth  has  several  undesirable  features; 

(1)  It  allows  the  transmission  of  excessive  noise  $ 

(2)  It  causes  unwanted  saturation  by  raising  the 
internal  signal  level;  along  with  this  it 
increases  the  requirements  for  components 
operating  at  high  power  levels,  since  output  peak 
power  is  generally  associated  with  the  higher 
frequencies! 

(3)  It  complicates  the  compensation,  since  the  high 
frequency  dynamics  of  the  components  cannot  be 
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neglected  in  determining  their  transfer  function^ 
and  the  compensation  must  effectively  offset  these 
additional  components 0 
This  reference  presents  an  analytical  method  for 
minimizing  bandwidth  by  transposing  the  problem  to  one  of 
minimizing  the  mean-square  value  of  a  transmitted  signal. 
The  method  is  good  for  any  configuration  and  for  both  tran- 
sient and  stochastic  signals*   It  uses  a  bandwidth  testo 
For  the  test,,  a  stationary  stochastic  noise  signal 
(characterized  by  that  to  which  the  system  is  expected  to  be 
subjected)  is  used  as  the  input  to  the  system..   The  result- 
ant output  is  passed  through  a  filter  whose  rms  output  is 
measured.  A  standard  system  is  defined  with  variable 
bandwidth.   It  is  fed  the  same  noise  signal  and  its  output 
passed  through  a  filter.   The  rms  value  of  the  control  system 
filtered  output  is  correlated  with  the  control  system 
bandwidth  by  comparison  with  the  standard  system  filtered 
output.   The  bandwidth  of  the  standard  system  is  adjusted 
until  the  two  rms  filtered  outputs  are  equals  by  definition^ 
the  two  bandwidths  are  then  equal.   Only  one  restriction 
exists?   The  standard  system,,  along  with  the  noise  source 
and  filter „  must  produce  a  monotonically  increasing  rms 
output  with  increasing  bandwidth?  otherwise ,  minimizing  the 
control  system  filtered  output  would  not  necessarily 
minimize  its  bandwidth.   By  this  scheme  the  problem  of 
minimizing  bandwidth  is  equivalent  to  determining  that 
overall  system  transfer  (or  weighting)  function  that 
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minimizes  the  mean  square  value  of  its  filtered  output 
subject  to  the  constraints  required  by  the  performance 
index  of  the   normal   input   signal. 

It  follows   that   additional   specifications   are    required 
for  establishing  the   noise    source 9    the    standard  system  and 
the   filter q      A  frequently  chosen  noise    source    is  white 
noise   because    it  has   components   at   all   frequencies   and   is 
easily  handled  analytically,.      The   filter  is   chosen  with 
regard  to   the    cutoff  characteristics  desired  in  the    systems 
the   filter  must  have    a  finite   mean-squared  output   since    this 
is   the    signal   to  be  minimized.      This   reference   recommends 
and  uses   a  pure   differentiator  of   (n-l)th  orderj,    requiring 
system  cutoff  at   the   rate    of  l/w11     The    standard   systems 
recommended  and  used  are   the   simple   binomial  filter  or  the 
Butterworth  filter   (see   p.220fl   reference   2)    chosen  so  that 
there    is  produced  a  finite  value   of  me  an- square   filtered 
output .      Since   the   problem  is  now  transposed  to  that  of 
minimizing  the   noise    transmitted  through  the    system,,    subject 
to   the   constraints   imposed  by  the   performance    index,,    another 
problem  in  variational   calculus   re  suits »      The    solution 
proceeds   as   indicated   in  the     paragraph     about      saturation; 
The   LaGrangian  multiplier  is   introduced  which  eventually 
leads   to   a  Wiener-Hopf  equation. 

(d)    Stability 

All  of  the    above    techniques  have    system  stability 
inherent   in  them  because   of  the    "spectrum  factorization" 
technique   required  to  solve   the  Wiener-Hopf  equation. 
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There   is   another  matter  to  considers      These    techniques  use 
equivalent  cascade   compensation,    requiring  manipulation  back 
to   the   feedback  loop   after  G„  (s)    has  been  determined*      The 
possibility  exists   that,    if  the   plant  elements   are   not 
minimum  phase    (zeros   in  RHP),    G-  G (s)    introduces   a  pole    to 
cancel   the   RHP  zero,    and   thus   Gc  (s)    is  unstable    in  its   own 
righto      To  prevent   this   the    above    techniques   require    that 
the   plant  elements  be    stable*      This   is  not   a  restriction,, 
since    all  that   is  necessary  is   to  provide    some   feedback 
loop   around  the   plant  elements   to  make    them  stable,    and  use 
the   equivalent   transfer   function  as   the   fixed  elements   for 
the   rest   of  the   design* 

(e)    Disturbances 

The   above    line   of  reasoning  applies   to  disturbances 
also?      The    fact   that   they  are   applied  at   a  place    other  than 
the    input  merely  requires   that   they  be   manipulated  out   to 
the    input o      This  modified  input   is   then  used  in  the   design* 
6*5      Graham  &  Lathrop  Methods12 

These   methods   are   not  Analytical  Design  Methods,,   but 
their  objective    is   the    same:      The    Synthesis   of   "Optimum" 
Transient   Response*      This   reference ,    first  of  all,    seems   to 
have    about   the   most  complete   expose    in  print   on  the   merits 
of   the   various  performance    indices*      Because    of  the   expose 
it   is   concluded   that  the   best  index  is   Integral-of-Time- 
Multiplied-Absolute-Error    (ITAE).      Since    this   function  is 
not   analytic 5    strictly  analytical  methods  cannot  be    applied* 

An  almost  exhaustive    study  was  made    on  the    following 
criteria,    as   applied  to  a  second  order  system  with  a   step 
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input    (this  was   a  zero-displacement  error  system  whose    trans- 
fer function  was   of  the    form     %  U)  =  = — V 


o 

a© 


Ti«f|ete)|<te  •  11  -fte*W^ 


Plots   of  the   various   criteria  versus   damping  ratio  resulted 
in  the    following? 

(1)  P^   and  P4  dictate    a.^f-0  for  minimum  error  and  were 
rejected  on  this   ground© 

(2)  P3   called  fovJ^0S  but  had  very  little    selectivity; 
but   it  can  be   handled  analytically  and  can  be 
mechanized  on  a  computer. 

(3)  P2   and  P5   called  for  jf  a=?»*5  °    both  are   easily 
mechanized  for  the   computer;    both  have   better 
selectivity  than  P3  with  P5   giving  much  sharper 
selectivity. 

(4)  Pg3   Pio*plls    called  for^/^r.4-,8    and  gave   good 
selectivity.      All  were   rejected  because    they  are 
too  difficult  to  handle   either  analytically  or  by 
computer* 
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^2^39^5  were   retained  and  applied  to   a  zero  velocity 
error   system  of  form 

T?      J S"+  2  f  5  ~l and  thane. 


applied  to  higher  order   systems.      It  was   found  that  P5(ITAE) 
was   the   only  criterion  which  retained  any  selectivity*      It 
was 5    the re fore s    concluded  that   the   ITAE   criterion  was  best* 
(It  is   interesting  to  note    that  one   of  the   advantages   of  ISE 
criteria  offered  by  Newton  is   the    fact   that   it   is  not  very 
selective •    thus   allowing  wide   latitude   practically  in 
adjusting  the    system.,) 

Next   in  the    study,    a   set   of  limitations   is  established? 
The    servo  must  be    a  "duplicator"?    namely,    it  will  be 
subjected  only  to   a   step   input   and  it  will  have   zero-displace- 
ment error   (meaning  the   closed  loop   transfer  function  has 
unity  in  the   numerator)*     With  these   restrictions   the   servo 
is  merely  a  low  pass   filter*      The    servo  could ,    therefore, 
be   represented  by  standard  forms   of  either  a  binomial  filter, 
or   a  Butterworth  filter,    or  a  "minimum  ITAE"   filter.      It  is 
a  simple    algebraic  matter   then,    to   force    the   coefficients   of 
the    servo  characteristic  equation  to  equal  the   coefficients 
of  the    standard  forms*      It  was   found   that   the   minimum  ITAE 
filter  response   combined  the   good  response   time   of  the 
Butterworth  filter  with  the    smaller  oscillation   (and  over- 
shoot)    of  the   binomial  filter*      However,    in  the   attempt   to 
extend  the  method  to  zero-velocity  and  zero-accele ration- 
error  systems,    step- function  responses   to   the   minimum  ITAE 

filter  were   no  better  than,    and   in  some   cases,   worse    than, 

94 


the  binomial  filter o 

The  use  of  standard  forms  does  not  involve  solving  for 

the  roots  of  equations ,  plots  or  graphical  construe tionsp 

integration,,  or  inverse  Laplace  transformations-   As  the 

reference  suggests 9    it  is  a  true  synthesis  method,,  in 

that  it  leads  directly  to  a  description  of  the  required 

system  in  terms  of  its  design  parameterso  There  are  a 

large  number  of  servos  which  fall  within  the  imposed  res trie- 

tionso 

6*4  Methods  of  Zaborszky  and  Diesel' 
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Looking  behind  all  the  indirectness  of  feedback  control 
system  design,,  it  develops  that  actual  system  design  is  based 
on  the  one  ultimate  measure  of  performance?   How  much  error 
is  there  at  the  times  when  the  system  output  is  utilized? 
Zaborszky  and  Diesel  attempt  a  mathematical  formulation  of 
such  a  measure  of  performance e   In  contrast  to  most  standard 
techniques  which  concentrate  on  isolated  phases  of  performance 
such  as  transient  or  steady  state,  this  measure  unites  all  of 
these,,  favoring  none9  in  a  single  concepto  This  measure  not 
only  concerns  itself  with  the  amount  of  the  error s   but  also 
the  times  when  the  output  is  being  utilized* 

This  measure  is  formulated  by  the  following  line  of 
reasonings 

(1)  The  specific  environment  and  function  of  each  control 
system  will  set  a  specific  penalty  valuation  on  errors  of  a 
given  sizeo   This  penalty  valuation  can  be  expressed  in  the 
form  of  a  penalty  function  P(e)s  a  single-valued  function  of 
the  error  e«  Error  itself  is  a  function  of  time,,  and  the 
valuation  may  vary  with  extraneous  parameters p  or  with  time*, 
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Therefor*  let  F<«)  -  F(<^+),  t,  TJT, .  "Ul  •  •  •  "Wl) 

(2)  A  second  element  of  performance  evaluation  is  the 

time  when  the  output  is  utilized*   For  simple  systems,*  all 
times  of  utilizing  the  output  are  generally  equally  probable, 
but  this  is  generally  not  so  in  advanced  or  complicated 
systems o   This  means  that  the  times  elapsing  from  activation 
of  the  system  to  all  times  of  utilization  of  its  output  can 
be  arranged  into  a  probability  distribution  p(t)» 

The  measure  can  be  mathematically  expressed  ass 


3   =  JoF(e)  fit)  dt 


where  the  symbol  \    is  the  end  sigma,  the  form  of  letter  sigma 
used  at  the  end  of  Greek  words „      For  a  deterministic  input, 
the  above  gives  the  average  value  of  the  penalized  error  at 
all  times  when  the  output  is  utilized*   If  the  environment 
consists  of  several  inputs «  or  is  stochastic ,  then 

where  the  bar  denotes  the  averaging  over  the  ensemble B   If 
p(t)  is  independent  of  the  input  (and  it  frequently  is), 
then 


5  = 


>o 


'o 

Conceptually  then,  the  above  integrals  are  the  average  value 
of  the  penalized  error  at  such  times  when  the  output  is 
utilized,  and  are  referred  to  as  the  "probabilistic  error" 
or  "end  sigma  error"*   As  a  performance  index,  the  integral 
is  called  the  "end  sigma  criterion" •   This  measure  unites  in 
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a  single  concept  the  transient  and  steady  states  of  operation 
as  well  as  any  intermediate  states*   None  of  these  states  is 
discriminated  for  or  against ,  because  p(t)  determines  the 
weight  allotted  to  each  state*   The  most  fascinating  thing 
about  the  above  integral  is  its  complete  generality?   All  of 
the  other  performance  indices  mentioned  in  this  paper  are 
simply  special  cases  of  this  integral* 

It  is  possible  to  evaluate  this  integral  in  the  S 
(frequency)  domain  if  the  functions  involved  are  Laplace 
transformable*   This  means  it  can  be  used  directly  with  a 
root  locus  plot  without  the  necessity  of  going  to  the  time 
domain*   It  can  also  be  evaluated  in  the  time  domain  in  the 
same  manner  as  indicated  in  Chapter  4*  Use  of  this  criterion 
for  optimizing  a  system  leads  again  to  the  Wiener-Hopf 
integral  equation*   The  second  reference  indicated  gives  a 
complete  design  process  to  obtain  the  impulse  response  function 
for  the  system  which  has  the  smallest  average  square  error 
for  such  times  when  its  output  is  used*   All  of  the  necessary 
equations  are  presented  in  an  organized  (albeit  complicated) 
form  for  the  programming  of  a  digital  computer*   As  they 
indicate p  the  solution  of  problems  of  a  complexity  common  in 
control  engineering  practice  today,,  requires  computers* 
6*5  Methods  of  Schultz  and  Ride out 

These  writers  propose  a  general  integral  of  the  form 


O 

which  is  the  same  as  that  of         Zaborszky  and  Diesel 
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with  no  consideration  to  the  times  of  utilization  of  the 
output.   However.,  their  attitude  is  the  recognition  of  the 
fact  that  physical  systems  can  never  respond  instantaneously., 
and  therefore  the  immediate  initial  error  should  be  ignored 
by  a  suitable  delay.   They  express  error  as  follows t 

e[±/x)  =  re  (*-*)-  c6t) 

where  r(t-tT)  represents  the  input  delayed  by  an  amount  °T  p 
and  c(t)  is  the  output.  This  error  is  called  "delay- error" 
and  it   is   the   one  used  in  the    indices   as  follows? 

Their  study  shows   that  the    integral  can  be  used  with 
transient   or   stochastic    inputs p    and   that  an  optimum  value 
of  V  exists   for  a  given  set  of  parameters.      They  have    also 
shown  that  IADE    (Integral  of  Absolute   value    of  Delay-Error) 9 
ISDE    (Integral   of  Square    of  Delay-Error) ,    and  ITADE    (Integral 
of  Time -weighted  Absolute   value   of  Delay-Error)    can  all  be 
easily  instrumented  on  an  analog  computer© 
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CHAPTER  7 
CONCLUSIONS 

7ol  General 

It  is   impossible   to  formulate   in  detail  a  universal 
approach  to  all  servo  design  problemsj,   but   it  is  possible 
to  list  in  proper   sequence   a  number  of  design  steps  which 
can  serve   as   a  guide •     As  mentioned  in  the   beginning., 
feedback  control   systems  have   penetrated  all  walks   of  life© 
All  kinds   of  engineers   and  mathematicians   are   in  the   fields 
all  developing  methods   to  fit   their  own  particular  problems • 
There    is   another  large   group  of  persons   in  this  field— 
this   group  never  really  designs   servos  per  se—they  look 
to   the   field  for  stimulation  and  challenge   to  their 
ingenuity?      the  mathematics   and  graphical  and  block  diagram 
manipulations   are   nfun"  <,    and  a  major  portion  of  them  are   not 
hard  to  learn.     And   so  long  as   the   design  is   a  paper  design^ 
no  "hardware M  experience   is  really  necessary*     As  automation 
makes   further  strides9    less  experience   is  needed,,   in  one 
sense j,   because  more   and  better  pieces   of  hardware   become 
available.     Every  day  it   is  more  nearly  possible    to  realize 
a  compensation  which  was  not  realizable    a  few  years   ago. 
But  at  the    same    time,,    two  other  things   are  happenings    1) 
cost  control  becomes  more   important ?    and  2)    advanced 
systems  used  with  refined  components  require   a  more 
sophisticated  analysis*) 

Given  a  set  of  specifications 9   many  solutions   to  the 
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design  problem  are  possible.   Owing  to  the  great  variety 
of  techniques  and  problems,  experience 9    in  the  final 
analysis 9   plays  just  as  important  a  role  as  the  use  of 
methodso   As  Gille   puts  it,  the  principle  most  often 
overlooked  is  that  a  feedback  control  system  constitutes 
an  entity,  and  each  component  of  it  must  be  considered  as  a 
part  of  the  overall  system,.  We  would  presume  that  this 
principle  never  really  strikes  home  until  one  has  attended 
the  school  of  whard  knocks". 
7.2  Analytical  Design  Methods 

There  can  be  no  doubt  that  these  techniques  are  rapidly 
coming  into  vogue.   The  general  reason  for  this  is  that 
higher  performance  is  continually  being  demanded?   Feedback 
Control  Systems  accomplish  all  kinds  of  sophisticated  jobs 
that  were  not  done  ten  years  ago  or  even  conceived  twenty 
years  ago.   The  rapid  development  of  improved  components 
continues  to  make  it  easier  for  the  designer  to  translate 
a  complicated  paper  design  into  a  useful  physical  system. 

The  first  hurdle  in  the  use  of  these  techniques  is 
to  bring  oneself  to  accept  the  idea  that  a  performance 
index  can  properly  express  the  specifications  of  the  system. 
Once  this  has  been  accepted,  application  of  the  methods 
becomes  much  more  palatable.   The  only  two  criteria  In 
extended  use  are  Mean  Square  Error  (MSE)  or  Integral  Square 
Error  (ISE)  and  Integral  Time  Weighted  Absolute  Error 
(ITAE).   Par  purely  analytical  techniques  only  MSE  can  be 
used. 
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Since  the  objective  of  the  design  is  the  minimization 
of  the  performance  index  j,  a  particular  advantage  of  these 
methods  lies  in  their  ability  to  recognize  inconsistent 
specifications^  in  the  sense  that  if  the  desired  performance 
index  is  less  than  the  theoretical  minimum  value 9   then  the 
specification  cannot  be  fulfilledo 

Another  advantage  lies  in  the  fact  that  the 
procedures  are  readily  susceptible  to  modification  to 
restrict  the  bandwidth  or  to  limit  saturation  tendencies o 

Another  advantage  lies  in  the  methods"  ability  to 
handle  virtually  any  kind  of  input  or  output.  Whereas 
most  Trial  and  Error  procedures  hinge  on  the  desired  output 
being  equal  to  the  input,,  these  methods  are  generally 
indifferent  as  to  what  the  desired  output  is©   This  means 
they  are  readily  able  to  accommodate  noisy  signals*   For 
that  matter^  they  can  handle  any  stationary  stochastic 
signal  as  long  as  it  is  possible  to  obtain  the  correlation 
function,,  or  correlation  function  transform© 

Another  advantage 9   albeit  a  philosophical  one,,  is  that 
the  mathematics  associated  with  the  methods  is  the  same 
type  of  mathematics  used  in  information  and  statistical 
theory  and  in  all  advanced  communication  theory©  Persons 
schooled  and  trained  to  think  in  those  lines  can  readily 
adapt  their  thinking  to  the  Analytical  Design  Methods* 

Besides  the  indicated  disadvantage  of  only  one 
performance  index0  another  disadvantage ©  that  could  be 
important  some  time s,  is  the  large  number  of  numerical 
calculations  associated  with  these  methods©   In  one  sense D 
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this  is  not  serious  sine©  any  method  requires  numerous 
calculations  if  the  system  is  an  advanced  ones  be  side sp 
computers  are  readily  available  today  to  do  much  of  this 
detail© 

Prom  the  viewpoint  of  minimizing  MSB,,  these  methods 
bring  home  very  clearly  three  theoretical  performance 
limitations  on  linear  systemss 

1)  Noise  and  disturbances  make  it  impossible  for  a 
system  to  establish  equality  between  the  desired  and  actual 
outputs  which  means  the  minimum  MSE  is  not  zero  under  these 
conditions j,  nor  as  low  as  it  could  be  in  the  absence  of 
the  noise  or  disturbance e 

2)  The  best  compensation  cannot  overcome  the  effect 
of  a  pure  time  delay  in  the  fixed  elements©  This  means 
that  a  feedback  control  system  cannot  predict  the  future 
value  of  a  signal  with  zero  errorp  since  it  inherently 
must  operate  on  present  or  past  information* 

3)  There  is  no  way  to  eliminate  completely  the  effect 
of  a  non°minimum=phase  fixed  element.   Intuitively  the 
cancellation  of  a  zero  in  the  RHP  can  be  approximated  but 
not  completely  accomplished  because  of  the  threat  of 
instabilityc   Thus  zero  MSE  cannot  be  obtained* 

We  know  that  a  practical  performance  limitation  is 
that  of  saturation,.   In  the  Trial  and  Error  Methods^ 
there  is  a  tendency  to  overlook  this  very  important  point 
until  one  attains  a  great  deal  of  experience©   In  the 
Analytical  Design  Methods*,  saturation  can  be  made  to  appear 
in  the  forefront  from  the  beginning  of  the  design© 
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7.5  Trial  and  Error  Design  Methods 

Categorically  the  transfer  function  approach  will  most 
directly  predict  frequency  response  while  the  root  locus 
approach  will  most  directly  predict  transient  responseo 
Use  of  the  system  equation  most  directly  predicts  stability, 
However 9    all  approaches  are  used  in  actuality* 

One  of  the  drawbacks  of  frequency  response  techniques 
is  the  difficulty  of  trying  to  visualize  the  transient 
response.   The  root  locus  method  fills  this  gap  because  the 
motion  of  the  closed-loop  poles  can  be  easily  observed  as 
the  gain  factor  is  varied.   The  root  locus  technique  then 
serves  as  a  fine  educational  tool  also*  Another  advantage 
lies  in  the  fact  that  the  compensation  employed  or  being 
investigated  is  readily  evident  and  the  problem  of 
unrealizable  compensation  seldom  presents  itself.  Direct 
determination  of  stability  is  also  easy  whereas  the 
frequency  response  techniques  can  sometimes  be  misleading 
for  indication  of  stability. 

A  peculiar  advantage  of  frequency  response  techniques 
(or  steady  state  analysis)  is  that  knowledge  of  the 
mathematical  model  is  not  a  requirement  and  it  is 
particularly  easy  to  select  a  compensator  to  cause  the 
response  curve  to  take  on  the  right  shape.   In  root  locus 
techniques  (and  for  that  matter5  in  Analytical  Design 
Methods)  j,  the  mathematical  models  of  different  parts  of 
the  servo  must  be  known^  in  practice  it  may  be  difficult 
to  determine  the  constants  required  for  such  a  representa- 
tion. 
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The  major  drawback  to  the  root  locus  is  that  accurate 
plotting  of  the  locus  is  a  time  consuming  task.  As 
indicated  earlier,,  a  number  of  theorems  for  approximating 
the  locus  are  available j  a  number  of  computing  schemes  have 
also  been  developed  including  one,,  by  one  of  the  authors,, 
for  use  on  the  NCR  102  compute ro 

Translation  from  frequency  response  to  transient 
response  is  also  difficult  at  best.  Here  again,,  digital 
computers  can  play  an  important  role  in  speeding  up  the 
arithmetic  required  in  the  various  approximation  methodso 

Analog  computers  of  course  greatly  facilitate 
analysis  and  design*   In  addition,,  good  analog  computers 
are  extremely  valuable  as  an  aid  to  remaining  in  the  linear 
zone,  or  observing  the  effects  caused  by  moving  out  of  the 
linear  zone.  To  recognize  their  use  in  this  respect  is 
merely  to  recognize  that  nature  just  is  not  linear* 
Probably  their  strongest  contribution  is  their  application 
to  solution  of  optimization  problems 9    after  the  configuration 
is  fixed. 

All  of  the  various  aids  and  charts  developed  in  these 
techniques  come  down  to  one  points   they  are  a  scheme  for 
relocating  the  roots  of  the  closed°loop  equation  by  jug- 
gling the  roots  of  the  open-loop  equation*   In  one  case 
the  given  specifications  define  a  closed-loop  equation  and 
the  problem  is  to  find  an  open-loop  equation  to  fit  it.   In 
the  second  more  advanced  case  the  problem  is  to  find  an 
open-loop  equation  in  which  the  plant  elements  fit  also. 
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7o4     Epilogue 

The    servo  designer  has   to  answer  three  major  questionss 

a)  Iffhat   is   the   Frequency  Response? 

b)  What   is   the   Transient  Response? 

c)  Is   it   stable? 

No  single  method  answers  these  questions*.  Careful  examina- 
tion of  a  proponent  of  a  particular  method  reveals  that  the 
method  works  well  for  some  particular  type  problem,,  but  has 
limitations  when  extended  to  some  other  problem© 

No  matter  what  method  is  consideredj,  the  first  step  is 
the  formulation  of  the  problem  by  gathering  the  appropriate 
specifications  and  plant  element  descriptions*.   Perhaps  the 
next  best  step  is  the  application  of  Analytical  Design 
Theory  to  the  development  of  an  appropriate  formula  for 
compensation,.   It  is  then  necessary  to  make  some  reasonable 
approximations  to  reduce  the  complexity  of  this  formula  or 
the  computations  associated  with  it,  after  which  the 
compensation  is  still  unduly  complex  and  generally  unre~ 
lizable*.  At  this  point  the  Trial  and  Error  Methods  should 
be  injected  into  the  amalgam  to  ameliorate  the  situation. 
Simpler  forms  of  compensation  found  by  these  methods  will 
almost  always  yield  performance  close  to  that  determined 
theoretically*  How  often  this  is  true  depends  upon  the 
broadness  of  the  minimum  of  the  chosen  performance  index*. 
Furthermore  several  different  forms  of  compensation  can  be 
found  and  compared  to  the  theoretical  one  for  their  relative 
efficiencies*,   The  reason  that  this  comparison  is  possible 
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is  that  it  is  relatively  easy  to  optimize  the  parameters 
after  the  selection  of  a  compensator  by  Analytical  Design 
Methods,,  since  the  system  becomes  one  of  a  fixed  configura- 
tion*  This  can  also  be  done  with  an  analog  computero 
However.*  the  practical  value  of  optimization  can  be  over- 
emphasized*  When  there  are  more  than  about  two  free 
parameters j,  optimization  is  "but  a  modern,,  more  systematic 
variation  of  the  very  old  engineering  practice  of 

-1  e 

compromise"  (Quote  from  Gille   )• 

Very  powerful  tools  are  available  today  for  synthesizing 
linear  feedback  control  systems..   All  are  of  value  9   providing 
the  designer  is  well  trained  in  their  use.   Apparently  no 
optimum  set  of  methods  exists,,  although  there  exists  a 
broad  enough  family  of  methods  to  arrive  at  "almost  the  best 
possible"  servo  for  a  given  set  of  conditions.   It  would  seem 
that  further  refinements  of  the  basic  methods  would  have  to 
await  additional  improvement  in  servo  components* 
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APPENDIX  A 
SELECTED  DEFINITIONS 

lo   Dynamic  Response  -  the  output  response  to  an  Input  that 

is  a  varying  function  of  time. 

20   Steady- St ate -Response  -  the  output  response  to  an  input 

that  is  constant  with  time* 

5e  Transient  Response  -  the  time  variation  of  one  or  more  of 

the  system  outputs  following  a  sudden  change  in  one  or  more 

of  the  system  inputs  or  the  derivatives  or  integrals  of  the 

system  inputs •  A  given  transient  response  must  be  referred 

to  the  type  of  input  that  caused  it. 

4.  Frequency  Response  -  the  variation  of  the  output  to  an 

input  which  is  a  const ant- amplitude  variable-frequency 

sinusoid* 

5e  Forced  Response  -  the  time  response  of  an  output  of  the 

system  to  an  arbitrary „   but  completely  defined,,  variation  of 

one  of  the  system  inputs.   Forced  response  is  distinguished 

from  transient  response  in  that  the  input  variation  associated 

with  the  forced  response  of  a  system  is  considered  as  a 

continuous  time  function  with  no  discontinuities  in  any  of 

its  derivatives©   (A  sinusoidal  input  is  a  special  case  of 

a  forcing  input  which  is  isolated  for  special  attention 

because  of  its  theoretical  importance ..) 

6.   Transient  Test  Inputs  - 

a.   Impulse  -  A  unit  impulse  is  a  time  function  that  is 
infinite  at  t^a  and  zero  everywhere  else.   It  is 
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defined  as  follows.,  where  ife(t°a)  is  a  unit  impulse 


function  occurring  at  t=a: 


i 

ri.(t-«)f(ojt=ftti) 


Step  -  the  unit  step  function  OLt  (t-a)  is  merely 
the  integral  of  the  unit  impulse  %Q   (t-a)o   It  is 
defined  as  follows^,  where  §l,(t°a)  is  a  unit  step 
occurring  at  tsa: 


e  -----    t 


o,  t<a 

Ramp  -   The  unit  ramp   function  ^,~(t-a)    is   the 
integral   of  the  unit   step  ^..(t-a).      The  unit  ramp 
is   defined  as   follows9   where  <§>     (t-a)    is   a  unit  ramp 
occurring  at   t=as  -£ 

S.A-k-a)   =f  £,(*-*)  dx 

d>  Parabolic  -  The  unit  parabolic  function  3„a(t-a) 
is  the  integral  of  the  unit  ramp  Sl  (t-a)«   ^iie 
unit  parabolic  is  defined  as  follows ,   where  c^.^ 
(t-a)  is  a  unit  parabolic  occurring  at  tsas 

J  o,  *  <;  a 

e«  Displaced  cosine  -  the  displaced  cosine  is  one  360° 
segment  of  a  cosine  wave  displaced  so  that  it  is 

always  positive  goings 
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Note  that  the  first  derivative  of  the  displaced 

2.TT 
cosine  is  zero  for  times  t-a  and  t-a-f*  — r-j 

Note  that  all  of  the  above  functions  are  equal  to 
zero  for  all  t(as  and  that  they  are  discontinuous^ 
or  one  or  more  of  their  derivatives  are  discontinu- 
ous at  the  instant  of  occurrence . 
7«>   The  Convolution  Integral 

a„   If  y(t)  is  the  inputs  x(t)  the  output^  and  w(t)  the 
impulse  response  of  the  system^,  then  the  output  x 
can  be  found  by  evaluating  the  convolution  integral 


or 


H&)s      wMJtfWjcrt, 


be   If  the  system  being  studied  is  a  physical  system, 
then 


U)(i)  =  0      "for  i  <  O 


and  the   convolution  integral  reduces   tos 
Do 


or  ° 

If  y(t)  and  w(t)  are  both  zero  for  t^Oj,  then  the 
convolution  integral  reduces  to 
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8,   The  Fourier  Transform  -  The  Fourier  transform  of  a 
function  and  its  inverse  are  defined  as  follows: 


where  o  =  the  complex  variable  cr+iuj 

The  Fourier  transform  is  applicable  to  functions  that  exist 

for  all  time  t«  To  insure  the  existence  of  the  Fourier 

transform  of  a  function,,  Diriehlefs  conditions  must  be 

satisfied* 

9e   The  Laplace  Transform  -  The  Laplace  transform  of  a 

function  and  its  inverse  are  defined  as  follows? 

where  s  s  the  complex  variable^-hiu^  •  Note  that  the  Laplace 
transform  is  used  for  functions  that  are  zero  for  t^O.   The 
constant  c  is  used  in  the  inverse  as  a  convergence  factor 
that  enables  one  to  apply  the  Laplace  transform  to  functions 
whose  Fourier  transforms  do  not  existe  A  function  must 
also  satisfy  the  Dirichlet  conditions  to  be  Laplace 
transformable  o 

10 ■   Gain  «  Gain  of  a  system  or  element  is  the  ratio  of 
magnitude  of  the  output  with  respect  to  the  magnitude  of 
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sinusoidal   input a      The    frequency  and  conditions   of  operation 

and  measurement  must  be    specifiedo 

11 o      Nyquist  Diagram  -   The   Nyquist  Diagram  is   a  closed  polar 

plot   of  a  loop   transfer  function  from  which  stability  may  be 

determined*     For  a  single-loop   system,,    it  is   a  map  on  the 

F(s)    plane   of  an  s-plane   contour  which  encloses   the  entire 

right  half  of  the   s-plane 9   excluding  poles   for  the   loop 

transfer  function  which  lie   on  the    imaginary  axis<> 

12 «>      Response   time   -   Response   time   is   the   time   required  for 

the   output  first  to  reach  a  specified  value   after  the 

application  of  a  step   input  or  disturbance • 

13 o      Rise   time   -   Rise   time   is   the    time   required  for  the   out= 

put   to  increase    from  one    specified  percentage   of  the   final 

value    to  anotherj,    following  the   application  of  a  step  inputo 

Usually  the    specified  percentages   are   10$  and   90$ B 

14 o      Settling  Time   -   The    settling  time   of  a  system  or 

element   is   the   time   required  for  the    absolute  value   of  the 

difference   between  the   output  and  its   final  value    to  become 

and  remain  less   than  a  specified  amounts    following  the 

application  of  a  step   input  or  disturbance  o      The    specified 

amount  is   often  expressed  in  terms   of  per  cent  of  the   final 

value o 

15 o   Gain  Margin  -  Gain  margin  is  the  amount  by  which  the 

magnitude  of  the  loop  ratio  of  a  stable  system  is  different 

from  unity  at  phase  crossover;  it  is  usually  expressed  in 

decibels<> 

16»  Phase  Crossover  -  Phase  crossover  is  a  point  on  the 
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o 
plot  of  loop  ratio  at  which  its  phase  angle  is  180  . 

17 .   Phase  Margin  -  Phase  margin  is  the  angle  by  which  the 

phase  of  the  loop  ratio  of  a  stable  system  differs  from 

180°. 

18 o   Gain  Crossover  -  Gain  crossover  is  a  point  in  the  plot 

of  loop  ratio  at  which  the  magnitude  of  the  loop  ratio  is 

unity* 

19,,   Loop  Ratio  -  Loop  Ratio  is  the  frequency  response  of 

the  primary  feedback  to  the  actuating  signal©  Under  linear 

conditions 9    the  ratio  is  expressed  as  GH  where  G  represents 

the  forward  elements  and  H  the  feedback  elements© 
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APPENDIX  B 
THE  FORMULATION  OF  THE  TRANSLATION 
AND  CORRELATION  FUNCTIONS 

B-l  Probability  Density  Functions 

The   analysis   of  stochastic   signals  requires   the  use   of 
probability  density  functions   and  other  statistical 
characterizations   such  as   the   average  value 9    the   root-mean- 
square  value   or  mean-square  value 9    and  the   correlation 
functions. 

The   probability  density  functions   are   direct  measures   of 

the   chance   of  occurrence   of  certain  events   in  a  process. 

The   first  probability  density  function  of  a  stochastic    (or 

> 

random)   variable  v(t)    is   denoted  by 

^1^v1p*i)  A  probability  that   the  variable   has   a 
value  v1  at  time    t. 
The    second  probability  density  function  is  denoted  by 

P2(vl9t1|V2t2)= probability  that   the   variable   has 
a  value  v^   at  time    t-^  and  a  value  v2  at   time   tg   simultaneously. 
For  a  stationary  stochastic   process    (one  whose    statistics 
are   independent  of  time),    ^(v^)    is   independent  of  time   t-jj 

p2(vls"tisv2pt2)    is   a  function  only  of  the    time   difference 
(t2-t^).      (Note   then,,    that  a  process  can  be   defined  with  one 
less  variable    if   it   is   stationary). 

In  general.,    the   average   or  mean  value   of  a  stochastic 
variable  v(t)    is  given  by    ^ 
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The   me  an- square   value  r  tX> 

is   given  bys  \/?(t)   ^        V°"  P(^J  ^  V 

The   root-me an- square    (rms)   value    is   given  by  the    square 
root  of  the  me  an- square   value*      The  variance   of  a  stochastic 
process   is   given  by 


V    £    [v-v] 


The  standard  deviation  0-  is  the  square  root  of  the  variance. 
It  can  be  expressed  as  follows? 


<r=  j>  -  ?ZJ 


Since  the  statistics  of  a  stationary  stochastic  variable 
are  independent  of  time,  the  mean  value  iss 

and  the  mean- square  value  is  given  by 

Two  commonly  used  probability  density  functions  are  the 
normal  distribution  and  the  Poisson  distribution*   The  normal 
distribution  is  given  bys 

where  P(v)dv  is  the  probability  of  finding  v  between  v  and 
v  +  dv. 

The  Poisson  distribution  is  given  by 
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n! 

where   P(  10.     At)    is   the   probability  of  finding  tt    events   in 
a  time   interval    At9    and  V    is   the   average   frequency  of 
occurrence   of  the   events • 


B-2     Correlation  Functions  for  Stationary  Stochastic   Signals 
The   autocorrelation  function    {fi.    C^Cj       ot  a  stationary 
stochastic  process  v(t)    is  defined  as   the  mean  value   of  the 
product  of  function  v  at  time   t  by  the   function  v   at  time 
(t  +*C  ): 

A  function  analogous   to   the   autocorrelation  function 
for  a  single    signal   is   the   cross-correlation  function  for  a 
pair  of   signals.      The   cross   correlation  function      yVU  C*v 
between  two   stationary  stochastic  processes  v(t)    and  u(t) 
is  defined  as   the  mean  value   of  the   product  of  the   function 
v  at  time    t  by  the   function  u  at  time    (t  Vf): 

Some   of  the  useful  properties   of  correlation  functions 
ares 

(1)    ^U  fr)    -  %v  (r°t)    (oven  function) 

(3)       jU*vo 


^tt)^ 


(p  or  a  stationary  signal  this  means  the  autocorrela- 
tion function  approaches  the  square  of  the  mean 
value  of  the  signal  as  *C  approaches  infinity.) 
(4)    Tv/u.  W  =  TUV  C"  V        (no,fc  an  oven  function) 

(5>    ll(t)K^'«)^W' 

(6)     ^«^u^=v"-IF 

We  observe  from  the  definition  of  the  autocorrelation 
function,,  that  the  mean  square  value  of  the  signal  equals 
the  value  of  the  corresponding  autocorrelation  function  with 

zero  argument:      ^f\>v  (O)  ~  ^^ 

B-3  Examples  of  Correlation  Functions 
(a)  Example  1 

A  common  type  stochastic  variable  used  is  the  case 
of  v(t)  as  a  rectangular  wave  with  values  4-  B   and  -M 
and  with  zero  crossings  located  at  event  points  that 
are  Poisson-distributed  in  time.   This  autocorrelation 
function  can  be  derived  ass 

-a* 


<GwCt)=p*6 


(b)  Example  2 

v(t)  is  a  rectangular  wave  with  amplitude  values 
distributed  in  any  fashion  and  with  zero  crossings 
Poisson-distributed  in  times 

where  (T  is  the  standard  deviation  of  the  amplitude 
distribution  and  V  is  the  mean  value  of  the  amplitude 
distribution.        ,,g 


(c)  Example  3 

v(t)  is  a  train  of  identical  finite  pulses  whose 
starting  points  are  Pois son- distributed  in  time  (known 
as  "shot  noise" )t  ^ 

where  f(t)  is  the  waveform  of  a  single  pulse  and 

•  DO 


V  =N  J   +-("tj  d"t.       (This  derivation  is   an  extension  of 
Campbell's  Theorem-   see   p«102s   reference   2) 
(d)   Example   4 

v(t)    is   pure   or  white   noises 

where  0    is   a  constant  that  depends   on  how  the   process 
is   generated^,    and    ^ft(t)    is   an  impulse.      For  example, 
if  white   noise    is   considered  as   a  limiting  case   of  shot 
noise   generated  by  exponential  pulses   of  amplitude   A9 
time   constant  Ts    and  area  s  under  the   pulse 9    then 


0   =    — 7—      ^5 


PC 

a 


B-4     Correlation  Function  Transforms 

Because   correlation  functions   are   completely 

defined  as   functions   of  a  time  variable   *C  9    they  are 

Fourier  transformable.      By  conventions,      ^=p  times   the 

Fourier  transform  of  a  correlation  function  is  called  a 

power  spectrum  or  a  power-density  spectrum*      Since    the 
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correlation  functions  involve  convolution  of  one  or  more 

functions;,  this  means  that  we  may  multiply  the  functions 

in  the  frequency  plane  and  this  is  the  significant  advantage 

of  taking  their  transforms. 

The  power  density  spectrum  3     ($)   of  a  stochastic  process 

is  defined  as 

The   cross-power- density  spectrum    ^     (s)  between  two 
stochastic   processes  v(t)    and  u(t)    is   defined  as 

The    inverse    transformations   are 

£>0 


Since  the  above  integrations  are  along  the  imaginary  axis 
of  the  s-planej,  the  inverse  transform  may  be  rewritten  in 
terms  of  the  real  frequency  U)  s 

*  no 
Recalling  that  the  mean- square  value  of  the  signal  is  the 
value  of  the  autocorrelation  function  with  ^T  a  O  9    then 


V-  DO 


DO 

Because  the  correlation  function  is  an  even  function  of  *C? 
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$w^U))can  bo   written  as  $vv(Uj)  »    or  ^(o)  =  I  Svv(u/J«/U, 
This  means   that   the    area  underneath  the   frequency  function 

5,  ^(Uj)      over   the    infinite   frequency  range    (commonly 
accepted  definition  of  power- density  spectrum)    is  equal   to 
the   mean-square  value   of  the   signal,. 

Useful  properties   of  the   power  spectra  are 

l>vv(2>)  -3>vV(-S>)       (even  function) 

B-5  Translation  Functions 

Since  transient  type  signals  exhibit  no  statistical 
properties ,  they  are  not  subject  to  description  by  correla- 
tion functions*   Newton  introduces  the  use  of  translation 
functions  (see  page  51*,  reference  2)  for  this  type  signal. 

If  x-^(t)  is  an  arbitrary  transient  signal 9    then  the 
autotranslation  function  X  y^j  ^s  defined  ass 

I„l-c)£  f  *.te)x.(++*)d± 

If  x^(t)    ans  X2(t)    are    two   arbitrary  transient   signals 9 
then   their  cross- translation  function    3I(aRJ  is   defined  as; 

xl2(«c)  £  |.]/.(-t)/»(-t+f)dt 

Although  these  functions  characterize  the  signals  this 
characterization  is  not  unique |  there  are  a  number  of 
different  functions  that  can  give  rise  to  the  same  trans- 
lation function0  Useful  properties  of  the  translation 

function  are  . 

X  (ftj  =  1,,  (-TJ   (even  function) 
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^X„io)^>llUi 


Note  from  the  definition  of  the  auto translation  function, 
that  if  t^O;  ,  t^ 

x, 

whore    1^   represents   the    integral- square   value    of  the   time 

function  x-^(t) 

B-6  Translation  Function  Transforms 

Because   of  Parse val's   Theorem   (see   p.  44,    reference   2)  9 
the    integral   1-^    stated  above,   can  be   expressed  in  terms   of 
its   transform  as 

The   Fourier  transform  of  a  translation  function    3T<avTj  is 
defined  as   follows? 


-ft© 


The  inverse  transform  is  defined  as 
Applying  this  definition  to  I  g 

X,   =    l„(0J  -    ^   J^po 

It  can  be  shown  (see  p0  57,  reference  2)  that 

X„(s)  =  Xi(-s)X,Cs), 

The  transforms  of  translation  functions  are  sometimes 
called  energy  density  spectra  as  contrasted  to  power  density 
spectra  in  correlation  functions«   It  is  to  be  noted  that  the 
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translation  functions  for  many  commonly  encountered  transient 
signals  are  infinite  since  their  defining  integrals  do  not 
converge.   Particular  examples  are  the  autotranslation 
functions  for  the  step  function  and  the  ramp.   Fourier 
transforms  for  these  functions  also  do  not  exist .   In  such 
situations  the  introduction  of  a  convergence  factor  will 
frequently  permit  a  solution  to  be  obtained.  However^  in 
general,,  the  determination  of  the  integral  square  value  of 
a  function  is  usually  done  in  the  time  domain  unless  the 
Fourier  transform  is  a  rational  function. 
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